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ABSTRACT 


The  objective  of  this  investigation  was  to  conduct  a  parametric  study  of 
the  effect  of  a  nonuniform  temperature  field  on  the  system  behavior  of  a 
trilayered  medium.  In  particular,  the  resulting  shear  and  normal  stresses 
along  the  media  interfaces  are  analyzed.  A  finite  element  model  utilizing  a 
recently  developed  element  which  provides  for  both  axial  and  lateral 
displacement  continuity  is  employed.  First,  the  effect  of  the  material 
properties,  that  is,  Young’s  Modulus  and  coefficient  of  thermal  expansion,  is 
examined.  Then,  the  effect  of  the  geometric  properties,  that  is,  length  and 
thickness  dimensions  of  the  midlayer,  is  analyzed.  Finally,  a  study  of  the 
effect  of  a  nonuniform  temperature  field  on  the  trimaterial  medium  is 
conducted. 
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I.  INTRODUCTION 


The  history  of  thermoelastic  stresses  on  layered  materials  begins  with 
Timoshenko.  In  a  1925  paper  [Ref.  1],  he  first  discussed  the  case  of  a  directly 
bonded  bi-material  configuration.  He  studied  bi-material  thermostats  modeled 
in  the  form  of  high  aspect  ratio  beams  subjected  to  uniform  heating  or  cooling. 
Bearing  stresses  in  closed  form  expressions  were  presented  and  it  was  noted 
that  there  was  not  an  elementary  way  of  determining  the  shear  stress 
distribution  along  the  bearing  surface.  The  curvature  of  the  beams  at  a  given 
temperature  was  shown  to  be  related  to  the  beam  geometry,  the  elastic 
constants  of  the  materials,  and  their  coefficients  of  thermal  expansion. 

Goland  and  Reissner  [Ref.  2]  and  Christensen  [Ref.  3]  studied  the  shear 
and  peeling  stresses  along  the  interface  of  single  lap  adhesive  joints  subjected 
to  axial  loads.  Goland  and  Reissner  derived  closed  form  analytical  solutions 
which  showed  that  the  peak  peeling  and  shear  stresses  occurred  at  the  joint 
edges.  Additionally,  it  was  shown  that  all  stresses  depended  on  the  joint 
dimensions  and  the  stiffness  properties  of  the  plates.  Christensen  performed 
his  analysis  using  the  finite  element  method.  He  determined  that  a  large 
stress  gradient  existed  across  the  thickness  of  the  adhesive  as  well  as  along 
the  lap,  especially  near  the  ends  of  the  lap.  Neither  study  investigated  the 
effect  of  temperature  on  system  behavior.  Burgreen  [Ref.  4]  analytically 
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computed  interfacial  shear  stresses  in  layered  beams  and  examined  the 
behavior  of  a  bi-metallic  beam  subjected  to  a  transverse  temperature  variation. 
Lau  [Ref.  5]  derived  closed  form  stress  analysis  solutions  for  a  chip  in  a  semi¬ 
infinite  substrate  and  a  chip  on  substrate  with  finite  thickness  subjected  to 
temperature  rises. 

A  typical  electronic  package  consists  of  a  semi-conductor  (chip)  attached 
to  a  substrate  material  by  an  adhesive  as  depicted  in  Figure  1.1.  Placing  more 
power  in  an  ever  decreasing  space  is  the  trend  in  electronic  component 
packaging.  Power  densities  of  12.5  MW/m3  have  been  reached  by  today's 
integrated  circuit  technology  [Ref.  6]  and  future  integrated  circuits  are 
expected  to  have  at  least  twice  this  power  density  leading  to  higher  operating 
temperatures.  The  resulting  high  operating  temperatures  give  rise  to  thermal 
cycling  causing  thermal  fatigue  at  the  different  material  interfaces  and  can 
lead  to  premature  failure  of  the  circuit.  To  prevent  failure,  a  maximum 
temperature  of  between  100°C  and  110°C  at  a  semiconductor  junction  has  been 
set  by  the  military  [Ref.  7]. 
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While  conducting  a  literature  search,  numerous  studies  concerning  the 
topic  of  thermal  response  of  electronic  packages  was  uncovered.  In  perhaps  one 
of  the  most  extensive  studies,  Suhir  [Ref.  8]  reviewed  work  conducted  on  bi¬ 
material  and  tri-material  electronic  packages  that  were  either  subjected  to  a 
uniform  temperature  or  an  applied  axial  load.  For  the  bi-material  configura¬ 
tion,  it  was  determined  that  the  maximum  peeling  and  shear  stresses  along 
the  interface  occurred  at  the  outer  edges  of  the  assembly.  Additionally,  it  was 
shown  that  the  stresses  depended  upon  the  axial  and  interfacial  compliances 
of  the  component  layers. 

Royce  [Ref.  9],  Riches  [Ref.  10],  and  Chung  [Ref.  11]  all  state  that  the 
primary  cause  of  mechanical  stresses  that  are  induced  in  electronic  packages 
are  a  result  of  the  difference  in  coefficients  of  thermal  expansion  of  the 
different  materials  in  the  assembly.  Royce  commented  that  thermal  cycling 
during  normal  operation  could  give  rise  to  thermal  fatigue  at  the  various 
material  interfaces  and  could  lead  to  premature  failure  of  the  structure.  Rich 
reported  that  the  types  of  stresses  that  arose  from  thermal  mismatch  consisted 
of  stresses  normal  to  the  chip  and  interfacial  stresses.  He  predicted  that 
normal  stresses  in  the  chip  were  at  a  maximum  in  the  chip  center  and 
decreased  to  zero  at  the  chip  edge.  However,  the  interfacial  stresses  were  zero 
at  the  center  and  rose  to  a  maximum  at  the  edge.  Rich  concluded  that 
compliance  and  thickness  of  the  chip  attachment  materials  were  critical  to 
controlling  the  stresses  imposed.  Two  ways  to  reduce  or  avoid  thermally 
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induced  stresses  in  electronic  packages  was  proposed  by  Chung.  The  first  was 
stress  reduction  by  coefficient  of  thermal  expansion  matching  and  the  second 
was  stress  reduction  by  the  use  of  extremely  flexible  adhesives. 

The  purpose  of  this  investigation  was  to  gain  further  insight  and  under¬ 
standing  into  the  effects  that  the  various  system  parameters  have  on  thermally 
induced  behavior  of  an  electronic  package  as  shown  on  Figure  1.1. 

To  achieve  this,  the  electronic  package  was  generalized  into  a  tri-material 
configuration  as  displayed  in  Figure  1.2  and  a  parametric  study  was  conducted 
using  finite  element  analysis.  First,  the  effect  of  the  material  properties, 
Young's  Modulus  and  coefficient  of  thermal  expansion,  was  examined.  Next, 
the  length  of  the  midlayer  was  varied.  Then,  the  effect  of  midlayer  thickness 
on  system  behavior  was  studied.  Finally,  nonuniform  temperature  fields  were 
imposed  on  the  tri-material  configuration. 


Figure  1.2  Tri-Material  Configuration 
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II.  FEM  MODEL  FORMULATION 


The  purpose  of  this  investigation  was  to  determine  the  stresses  resulting 
from  a  nonuniform  temperature  field  on  a  tri-material  configuration.  To 
accomplish  this,  a  finite  element  method  (FEM)  program  was  constructed. 
This  program  is  presented  in  Appendix  A  and  the  program  exec  is  in 
Appendix  B. 

A.  PROGRAM  DEVELOPMENT 

A  recently  developed  element  which  provides  for  both  axial  and  lateral 
displacement  continuity  was  utilized  in  the  stress  formulation.  As  shown  in 
Figure  2. 1,  each  element  has  six  degrees  of  freedom;  axial  displacements  at  the 
four  comer  points;  and  lateral  displacements  at  the  two  ends. 


y 


Figure  2.1  A  Typical  Element  With  Six  Degrees  Of  Freedom 


The  lateral  displacement  field  is  determined  by  linear  interpolation 

(2-1) 
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where  vt  and  v2  are  the  lateral  displacements  at  the  ends  and  where  the  linear 
shape  functions  Nx  and  N2  are  given  by 

(2.2) 

resulting  in 

vW=E^)v  1  (2.3) 

i- 1 

The  axial  displacement  field  is  assumed  to  be  linear  in  both  the  a  vial  and 
transverse  direction.  That  is, 

u*=Hfu^H[u[  (2.4) 

i (2-5) 
where  superscripts  b  and  t  refer  to  the  bottom  and  top  respectfully.  The  linear 
shape  functions  Hx  and  H2  are  given  by 

Hl-Ul-Z  (2.6) 

ft  n 

resulting  in 

2 

u(j:j)=£  W^)(//,0')U,‘.ffJ(j)«,‘]  (2.7) 

(-1 

From  strain-displacement  relations 

/  j  2  \ 

zx=^r  -He  N/WffjM#,'  <2  8> 

ax  (.i  ) 
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or 


and 


(2.9) 


du  dv 


(2.10) 


(2.11) 


or 


vE  Ni 


i-i 


M,  »  3H2  , 


t  + 


E§*i 

l-i  «* 


The  derivatives  of  the  linear  shape  functions  are 

dNt  i  dN2_i  dHx_  i  dH2  i 
dx  l  d x  l  dy  h  dy  h 


This  results  in 


(2.12) 


(2.13) 


(2.14) 


and 


v2-vi 


/ 


(2.15) 
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Defining  the  axial  displacement  vector  as 


W  u/  (2.16) 

and  the  stiffness  matrix  as 

0TJ  =  fJfoElBFdydx  (2. 17) 


where  the  B  vector  is 


{bF^h, 
dx  1 


dN . 

-s* 


aw,  aw,  . 

-s* 


and  the  force  vector  due  to  temperature  as 

gives  the  bending  matrix  equation  as 


(2.18) 


(2.19) 


(2.20) 


Equation  2.20  defines  the  bending  behavior.  Additionally,  the  stiffiiess  matrix 
can  be  shown  as 


2  1-2-1 

1  2-1-2 

-2-12  1 
1-212 


(2.21) 


Behavior  due  to  shear  is  obtained  as  follows.  Define  the  row  vector  of 
displacement  degrees  of  freedom  as 


{fi/^4*  u{  Vj  ih  Vj> 


(2.22) 
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The  shear  stiffness  matrix  is  given  by 


[jy  =  f‘fQk{B/\G(B')(fydx  (2.23) 


where 


.  A  .  dH.  dH- 

W4N'i£  N 


dN,  dH,  a H.  dN,, 

dr  dy  dy  dx 


(2.24) 


which  gives  the  equations  for  shear  behavior 


[jyfojMO} 

Additionally,  the  stiffness  matrix  can  be  shown  as 


/ 

-l 

1 

l 

-/ 

_1 

Ah 

Ah 

2 

Ah 

4A 

2 

-l 

l 

-1 

-l 

/ 

1 

Ah 

Ah 

2 

Ah 

4A 

2 

1 

-1 

h 

1 

-1 

-h 

2 

2 

l 

2 

2 

l 

/ 

-/ 

1 

/ 

-l 

-1 

Ah 

Ah 

2 

4/t 

Ah 

2 

-l 

l 

-1 

-/ 

l 

1 

Ah 

Ah 

2 

AH 

Ah 

2 

-1 

1 

-h 

-1 

1 

A 

2 

2 

l 

2 

2 

/. 

The  system  matrix  for  bending  and  shear  is 


(2.25) 


(2.26) 


(2.27) 
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Let 


a=2l  b=. a  c=M  d=° 

Ah  l  6 1  2 


(2.28) 


then 


a+2c  -a+c  d  a -2c  -a-c  -d 

-a+c  a+2c  -d  -a-c  a -2c  d 

d  -d  b  d  -d  -b 

a -2c  -a-c  d  a+2c  -a+c  -d 

-a-c  a -2c  -d  -a+c  a+2c  d 

-d  d  -b  -d  d  b 


(2.29) 


The  force  vector  due  to  a  temperature  field  is 


where 


{Fe}= Cl*f^f*{B}dydx 


(2.30) 


Q*=(£aA7y 


(2.31) 


JoJc 


oJo  aiv. 


(2.32) 
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Let 


(2.33) 


* 

then 


{F'}=A'< 


-1 

>1 

0 

1 

1 

0 


After  solving 


(2.34) 


(2.35) 


for  {8},  the  stresses  can  be  calculated.  For  the  e*  element,  the  bending  stress 
is 


or 


o  -Et 

X  X 


(2.36) 


(2.37) 


(2.38) 


where  aj™  and  a*  are  the  bending  stresses  acting  on  the  element  at  the 
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bottom  and  top  respectfully.  The  shear  stress  is 


(2.39) 


or 
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(2.41) 


where  xj™  and  are  the  shear  stresses  acting  on  the  e*  element  at  the 
bottom  and  top  respectfully.  Finally 


,  Gj(82+53Mfil  + 

'^=7f» — 


*,)  2(8^83) 


z 


(2.42) 


B.  PROGRAM  INPUT 

The  input  file  for  the  finite  element  program  is  displayed  in  Appendix  C. 
A  mesh  grid  with  200  elements  as  shown  in  Figure  2.2  is  generated.  The  input 
file  allows  for  variation  in  both  the  element  layer  thicknesses  and  the  element 
row  widths.  Values  for  Young's  Modulus,  Poisson's  Ratio,  and  the  coefficient 
of  thermal  expansion  can  be  inputted  row  by  row.  Additionally,  the 
temperature  for  each  level  of  the  mesh  can  be  inputted.  Appendix  C  is 
annotated,  showing  the  placement  of  the  different  inputs.  All  dimensions  are 
in  millimeters. 
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C.  PROGRAM  OUTPUT 


The  output  file  from  the  finite  element  program  is  shown  in  Appendix  D. 
The  force  and  displacement  vectors  for  each  node  are  outputted.  The  top  and 
bottom  bending  stresses  and  shear  stress  for  each  element  are  also  displayed. 
Additionally,  the  normal  stresses  along  the  material  A  and  material  B 
interface  and  the  material  B  and  material  C  interface  are  outputted. 


]VL 

vr 

E 

u 

E 

u 

I 

M 

U 

,  A 

i 

1 

1 

1 

1 

_ 

■ 

■ 

■ 

- 

Wk 

i 

■ 

'i 

i 

1 

i 

i 

1 

1AT 

ERIA 

L  B 

i 

1 

■ 

i 

1 

■ 

i 

1 

— 

i 

■ 

i 

■ 

i 

i 

i 

■ 

■ 

i 

l 

1 

IS 

1 

1 

■ 

i 

9 

1 

1 

1 

1 

1 

Figure  2.2  Finite  Element  Analysis  Mesh  For  Tri-material  Configuration 
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III.  EFFECT  OF  MATERIAL  PROPERTIES 


This  chapter  studies  the  effect  of  Young's  Modulus  and  the  coefficient  of 
thermal  expansion  on  the  tri-material  configuration  of  Figure  3.1.  Several 
computer  runs  were  made  assigning  different  Young's  Modulus  values  to  the 
three  materials  while  maintaining  a  constant  coefficient  of  thermal  expansion 
value.  No  stresses  resulted.  As  a  consequence,  the  remainder  of  this  chapter 
is  concerned  with  variations  in  the  coefficient  of  thermal  expansion  (a^.  Table 
3.1  shows  the  parameters  used  for  the  11  cases  that  were  evaluated.  For  all 
cases,  Young's  Modulus  was  given  a  constant  value  of  100  GPa  for  the  three 
materials  and  a  change  in  temperature  of  100°C  was  imposed. 


a - 10 - *| 


|  ALL  DIMENSIONS  IN  MILLIMETERS 

CENTERLINE 

Figure  3.1  Tri-material  Configuration 
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TABLE  3.1  COEFFICIENT  OF  THERMAL  EXPANSION  STUDY 


The  eleven  cases  were  grouped  into  four  combinations  to  allow 
comparisons  to  be  made.  The  combinations  were: 

•  Cases  1  and  2.  aA  and  ac  were  set  equal  and  aB  was  allowed  to  range 
from  one-half  a  a,  otc  to  twice  occ. 

•  Cases  3,  4,  and  5.  a A,  ctg,  and  were  given  different  values  in  different 
combinations. 

•  Cases  3,  6,  7,  and  8.  aB  and  ac  were  held  constant  but  at  different  values 
while  aA  was  allowed  to  range  at  values  greater  than  aB  and  otc. 

•  Cases  7,  9,  10,  and  11.  aA  and  aB  were  held  constant  but  at  different 
values  while  etc  was  allowed  to  range  at  values  less  than  aA  andaB. 
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A.  CASES  1  AND  2 


Cases  1  and  2  investigated  how  the  tri-material  configuration  responded 
when  the  top  and  bottom  layers  (materials  A  and  C)  held  the  same  value  of  the 
coefficient  of  thermal  expansion  while  the  middle  layer  was  given  values  lower 
(case  1)  and  higher  (case  2). 

Figure  3.2  is  a  plot  of  the  normal  stress  along  the  upper  interface  (between 
materials  A  and  B)  and  lower  interface  (between  materials  B  and  C)  as  a 
function  of  position.  For  case  1,  both  the  upper  and  lower  interfaces  are  in 
tension,  showing  that  -he  top  and  bottom  layers  with  their  larger  amount  of 
expansion,  wanting  to  peel  away  from  the  middle  layer.  In  case  2,  the  normal 
stresses  at  the  upper  and  lower  interfaces  are  all  compressive,  that  is,  have 
bearing  stresses.  The  middle  layer  wants  to  expand  more  than  the  top  and 
bottom  layers  and  pushes  against  them.  For  both  cases,  the  maximum  stresses 
occur  at  the  ends  of  the  configuration. 

Figure  3.3  is  a  graph  of  the  shear  stresses  acting  on  the  upper  and  lower 
interfaces.  The  shear  stresses  for  case  1  tend  to  axially  extend  the  midlayer 
while  the  shear  stresses  for  case  2  tend  to  axially  compress  the  midlayer. 
Figure  3.4  shows  the  shear  and  normal  stresses  acting  on  the  midlayer  for 
cases  1  and  2. 
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Case  aA(xlO  /°C)  c^UfTC)  a^xMTTC) 

Case  1  100  ~  50  100 

Case  2  100  200  100 

Figure  3.2  Normal  Stress  At  Interfeces  For  Cases  1  And  2 
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Figure  3.3  Shear  Stress  At  Interfeces  For  Cases  1  And  2 


B.  CASES  3,  4,  AND  5 

Cases  3,  4,  and  5  investigated  the  response  of  the  tri-material  configura¬ 
tion  when  all  three  layers  had  different  values  for  the  coefficients  of  thermal 
expansion.  Case  3  was  arranged  with  magnitudes  as:  aA>  aB>  ac,  case  4  as: 
otA>  ac>  aB,  and  case  5  as:  ocB>  aA>  ac. 

Figure  3.5  is  a  plot  of  the  normal  stresses  as  a  function  of  position  along 
the  upper  and  lower  interfaces.  Case  3,  with  the  coefficients  of  thermal 
expansion  in  decreasing  order  from  top  to  bottom,  displayed  minimal  peeling 
stresses  along  the  interfaces.  Case  4,  with  the  coefficient  of  thermal  expansion 
lowest  in  the  midlayer,  exhibited  characteristics  similar  to  case  1.  The  normal 
stresses  were  peeling  on  both  interfaces.  Case  5,  with  the  coefficient  of 
thermal  expansion  in  the  middle  layer  the  largest,  behaved  like  case  2.  The 
normal  stresses  along  the  upper  and  lower  interfaces  were  bearing  throughout. 

The  plots  for  the  shear  stresses  along  the  upper  and  lower  interfaces  are 
shown  in  Figure  3.6.  In  case  3  the  shearing  stresses  were  the  lowest,  with 
material  B  wanting  to  elongate  axially  along  the  upper  interface.  The 
midlayer  wanted  to  axially  elongate  for  case  4  and  to  axially  compress  for  case 
5.  This  was  the  same  behavior  that  was  observed  for  cases  1  and  2,  respectful¬ 
ly.  Figure  3.7  displays  the  shear  and  normal  stresses  acting  along  the 
midlayer  interface. 
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C.  CASES  3,  6,  7,  AND  8 

This  study  compared  the  responses  of  the  tri-material  configuration  when 
aB  and  occ  were  held  constant  at  values  lower  than  aA  (with  aB>  a, g),  and  aA  was 
then  allowed  to  increase  as  follows:  aA  (case  6)  <  aA  (case  7)  <  aA  (case  3)  <  aA 
(case  8). 

Figure  3.8  is  a  graph  of  the  normal  stresses  along  the  upper  and  lower 
interfaces  for  the  midlayer.  When  the  value  of  the  coefficient  of  thermal 
expansion  for  the  upper  layer  was  kept  reasonably  close  to  the  coefficient  of 
thermal  expansion  of  the  midlayer  (within  a  factor  of  two)  as  shown  by  cases 
6  and  7,  the  normal  stresses  tended  to  be  minimal  and  were  mainly  bearing. 
However,  when  aA  was  further  increased  as  in  cases  3  and  8,  the  magnitudes 
of  the  normal  stresses  increased  and  became  peeling  throughout  the  upper  and 
lower  midlayer  interfaces. 

The  shear  stress  plot  for  the  upper  and  lower  interfaces  is  given  in  Figure 
3.9.  Cases  6  and  7  resulted  in  shearing  stresses  at  both  ends  of  the  upper  and 
lower  interfaces  with  minimal  shear  stresses  along  the  majority  of  the 
interface  lengths.  In  cases  3  and  8,  the  shearing  stresses  on  both  interfaces 
tended  to  axially  stretch  the  midlayer  with  case  8  having  a  larger  magnitude 
of  stress.  Figures  3.10  and  3.11  display  the  shear  and  normal  stresses  acting 
on  the  midlayer  for  cases  3,  6,  7,  and  8. 
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A  =  aUUXlU  /  U _  V 379  MPa  aA=4UUXl 

ag=  lOOxlO^C  ac=  lOxlO^C 
Ea  =  Eb=Ec=  100  GPa 

Figure  3. 10  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  3  And  8 
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D.  CASES  7,9,  10,  AND  11 

In  this  grouping,  the  coefficients  of  thermal  expansion  for  materials  A  and 
B  were  held  constant  at  values  greater  than  material  C  (with  aA  >  otg),  and  otc 
was  allowed  to  increase  as  follows:  ac  (case  9)  <  ac  (case  7)  <  ac  (case  10)  < 
(case  11). 

Figure  3.12  is  a  plot  of  the  upper  and  lower  interface  normal  stresses  for 
the  four  cases.  Similar  characteristics  were  displayed  by  cases  7  and  9.  Along 
the  upper  interface  the  normal  stress  was  initially  peeling,  decreased  to  zero, 
and  then  became  bearing.  At  the  lower  interface  the  normal  stresses  were 
initially  bearing,  went  to  zero,  and  then  became  peeling.  For  cases  10  and  11 
the  normal  stresses  were  peeling  on  both  the  upper  and  lower  interfaces. 

A  plot  of  the  shear  stresses  acting  on  the  midlayer  is  shown  in  Figure 
3. 13.  As  can  be  seen,  along  the  upper  interface  all  four  cases  displayed  similar 
shear  stress  characteristics  with  the  surface  being  pulled  in  axial  tension. 
However,  at  the  lower  interface,  cases  7  and  9  showed  axial  compression  while 
cases  10  and  11  displayed  axial  tension.  The  shear  and  normal  stresses  acting 
on  the  midlayer  for  these  four  cases  are  shown  in  Figure  3.14  and  3.15. 

E.  SUMMARY 

A  review  of  the  peeling  and  shear  stresses  acting  along  the  upper  and 
lower  interfaces  reveals  force  equilibrium  in  both  the  lateral  and  axial 
directions.  However,  moment  equilibrium  is  not  satisfied  exactly.  This 
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Figure  3.13  Shear  Stress  At  Interfeces  For  Cases  7, 9,  10,  And  11 


deviation  is  most  probably  a  result  of  the  type  of  finite  element  that  was 
chosen. 

As  expected,  when  materials  A  and  C  were  given  the  same  properties, 
symmetrical  stress  conditions  resulted.  This  was  clearly  shown  in  cases  1 
and  2.  When  the  coefficient  of  thermal  expansion  in  the  midlayer  was  smaller, 
the  upper  and  lower  layers  wanted  to  separate  from  the  middle  layer.  When 
the  coefficient  of  thermal  expansion  for  the  midlayer  was  greater,  the  top  and 
bottom  layers  tended  to  compress  it. 

Maximum  peeling  stresses  occurred  when  the  coefficient  of  thermal 
expansion  for  the  top  or  bottom  layer  was  much  greater  than  the  middle  layer 
as  shown  in  case  4.  Similarly,  the  maximum  bearing  stresses  arose  when  the 
midlayer  coefficient  of  thermal  expansion  was  greater  than  the  upper  and 
lower  layers  as  evidenced  by  case  5. 

The  relative  difference  in  the  coefficient  of  thermal  expansion  between  the 
layers  played  a  similar  role  for  the  shear  stresses.  The  maximum  tendency  to 
axially  elongate  the  midlayer  occurred  when  the  coefficient  of  thermal 
expansion  for  the  middle  layer  was  much  less  than  the  upper  or  lower  layers. 
Additionally,  the  maximum  tendency  to  axially  compress  the  midlayer 
happened  when  the  coefficient  of  thermal  expansion  for  the  midlayer  was  much 
greater  than  the  top  or  bottom  layers. 
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Minimal  normal  stresses  resulted  when  the  values  of  the  coefficients  of 
thermal  expansion  for  the  upper,  middle,  and  lower  approached  each  other. 
The  conditions  for  minimum  shearing  stresses  were  the  same. 

The  study  of  cases  7,  9,  10,  and  11  produced  an  interesting  result.  As  can 
be  seen,  cases  7  and  9  exhibited  very  similar  stress  states  even  though  the 
coefficients  of  thermal  expansion  for  the  bottom  layer  differed  by  a  factor  of 
ten.  More  computer  runs  were  conducted,  further  reducing  the  coefficient  of 
expansion  by  several  factors  of  ten.  The  changes  in  the  state  of  stress  was 
minimal. 

Finally,  when  the  coefficient  of  thermal  expansion  for  the  tri-material 
configuration  was  kept  constant  and  uniform,  changing  of  the  Young’s  Modulus 
values  did  not  produce  any  resulting  thermally  induced  stresses. 
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IV.  EFFECT  OF  MIDLAYER  THICKNESS 


This  chapter  investigates  the  effect  of  changing  the  thickness  of  the 
midlayer  (material  B)  on  the  tri-material  configuration.  Five  different 
midlayer  thickness  cases  were  studied  and  compared.  As  displayed  in  Figure 
4.1,  the  following  thicknesses  for  material  B  were  used:  1.00mm  (case  5), 
0.50mm  (case  12),  0.25mm  (case  13),  0. 10mm  (case  14),  and  0.05mm  (case  15). 
For  all  cases,  the  dimensions  of  materials  A  and  C  were  kept  constant  (10mm 
length  and  1.0mm  thickness).  Additionally,  the  length  of  the  midlayer  was 
held  at  5mm  and  the  tri-material  configuration  was  given  a  temperature 
increase  of  100°C.  The  material  properties  for  the  five  cases  were  given  the 
following  constant  values:  EA=  EB=  Ec  =  100  GPa,  aA=  lOOxlO'V’C,  aB= 
SOOxlO'V’C,  and  ac  =  lOxlO^C. 

A.  NORMAL  STRESS  DISTRIBUTION 

Plots  of  the  normal  stresses  along  the  upper  and  lower  interfaces  of  the 
midlayer  for  the  five  cases  are  displayed  in  Figure  4.2.  Both  the  upper  and 
lower  normal  stress  distributions  exhibited  similar  characteristics.  The 
stresses  were  always  bearing  along  the  interfaces  for  all  five  cases  with  the 
maximum  stresses  occurring  at  the  left  and  right  ends  of  the  midlayer.  For 
every  case,  the  minimum  bearing  stress  was  developed  midspan  of  the 
interface.  As  the  thickness  of  material  B  was  decreased,  the  bearing  stresses 
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Normal  Stresses  At  The  Interfeces  For  Cases 


similarly  decreased.  However,  the  shape  of  the  normal  stress  distribution 
curves  remained  comparable. 

Comparisons  were  made  between  the  maximum  of  the  normal  stresses  of 
all  five  cases  on  the  upper  and  lower  interfaces  and  the  remainder  of  the 
normal  stresses.  Figures  4.3  and  4.4  are  two-dimensional  and  three-dimen¬ 
sional  plots  of  the  normal  stress/maximum  normal  stress  for  the  five  cases 
across  the  upper  and  lower  interfaces,  respectfully.  For  both  interfaces,  the 
maximum  normal  stress  occurred  in  case  5  at  the  right  hand  side  of  the 
midlayer.  Similar  characteristics  were  displayed  by  all  cases  for  both  the 
upper  and  lower.  As  the  thickness  of  the  midlayer  decreased,  the  ratio  of  the 
normal  stress  to  the  maximum  normal  stress  decreased.  Of  interest  was  that 
the  thinner  material  B  became,  the  extent  of  the  midspan  that  was  at  a 
minimum  grew.  Figures  4.5  and  4.6  show  the  bearing  stresses  acting  along  the 
upper  and  lower  interfaces  of  material  B. 

B.  SHEAR  STRESS  DISTRIBUTION 

Figure  4.7  contains  plots  of  the  shear  stresses  across  the  upper  and  lower 
interfaces  of  material  B  for  cases  5,  12,  13,  14,  and  15.  Cases  5  and  12  (with 
their  greater  thicknesses  of  material  B)  displayed  similar  attributes.  The 
shear  stresses,  upper  and  lower,  tended  to  axially  compress  the  midlayer. 
Cases  13,  14,  and  15  (with  midlayer  thicknesses  of  0.25mm,  0.10mm,  and 
0.05mm)  showed  peak  shear  stresses  at  the  left  and  right  ends  of  the 
interfaces,  but  along  the  majority  of  the  length  of  span  the  shear  stresses 
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1  Material  B  Height  (mm)  111  0,5  0.25  0.1  1  0.05 

aA=  lOOxlO^C  Og=  300x10‘6/°C  ac  =  lOxlO^C  EA  =  EB=EC=  lOOGPa 
Figure  4.5  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  5, 12,  And  13 
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_ Case  5  Case  12  Case  13  Case  14  Case  15 

Material  B  Height  (mm)  1  0.5 _ 0.25 _ 0J. _ 0.05 

aA=  lOOxlO^C  aB=300xl0‘6/°C  ac  =  lOxlO^C  EA  =  EB  =  EC=  lOOGPa 
Figure  4.6  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  14  And  15 
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fluctuated  around  zero.  This  fluctuation  of  positive  and  negative  values  of  the 
shear  stresses  from  point  to  point  may  have  been  due  to  a  breakdown  in  the 
finite  element  program  in  dealing  with  the  very  small  numerical  values 
present. 

Comparisons  of  the  ratio  of  the  shear  stress  at  a  point  to  the  maximum 
shear  stress  found  in  the  group  were  conducted.  Figures  4.8  and  4.9  are  two- 
dimensional  and  three-dimensional  plots  of  the  shear  stress/maximum  shear 
stress  for  the  five  cases  along  the  upper  and  lower  interfaces,  respectfully.  For 
both  the  upper  and  lower,  the  maximum  shear  stress  was  located  at  the  left 
end  of  material  B.  The  two  figures  were  similar,  with  the  shear  stress  ratios 
being  maximum  at  the  extreme  ends  of  the  midlayer  and  then  decreasing  to 
a  minimum  towards  midspan.  Of  interest,  like  the  normal  stress,  was  that  as 
the  thickness  decreased,  the  range  of  the  midspan  that  was  at  a  minimum 
increased.  Figures  4.5  and  4.6  display  the  shear  stresses  acting  along  the 
upper  and  lower  interfaces  for  the  midlayer. 

C.  SUMMARY 

The  effect  of  decreasing  midlayer  thickness  on  the  tri-material  configura¬ 
tion  was  displayed  in  a  number  of  ways.  The  normal  stresses  remained 
bearing  throughout  as  the  midlayer  became  thinner.  The  maximum  normal 
and  shear  stresses  decreased  with  decreasing  thickness.  Minimal  normal  and 
shear  stresses  extended  out  further  along  the  span  and  the  stress  distributions 
became  more  local  in  nature  with  the  thinning  of  material  B.  A  conclusion 
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Figure  4.8  2-D  And  3-D  Plots  Of  Shear  Stress/Maximum  Shear  Stress  At  The 
Upper  Interface  For  Cases  5, 12, 13, 14,  And  15 


46 


that  could  be  drawn  would  be  that  as  the  thickness  of  the  midlayer  continues 
to  decrease,  that  the  only  stresses  of  any  consequence  along  the  span  exist  at 
the  extreme  left  and  right  ends  of  the  interfaces.  The  normal  stresses  were 
only  one-third  to  one-half  the  shear  stresses.  For  both  the  upper  and  lower 
interfaces  the  maximum  normal  stresses  remained  at  the  far  left  edge  in  all 
cases.  The  maximum  shear  stresses  along  the  lower  interface  moved  from  the 
left  edge  to  the  right  edge  as  the  midlayer  thickness  was  decreased.  Finally, 
the  maximum  shear  stresses  along  the  upper  interface  remained  at  the  left 
edge  for  all  cases. 
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V.  EFFECT  OF  MIDLAYER  LENGTH 


In  this  chapter,  the  effect  of  the  length  of  the  midlayer  (material  B)  on  the 
tri-material  configuration  was  investigated.  Six  different  midlayer  lengths 
were  compared  and  studied.  As  shown  in  Figure  5. 1,  the  following  lengths  for 
material  B  were  used:  9mm  (case  16),  7.5mm  (case  17),  5mm  (case  12),  2.5mm 
(case  18),  1mm  (case  19),  and  0.5mm  (case  20).  For  all  cases,  the  dimensions 
of  materials  A  and  C  were  kept  constant  (10mm  length  and  1mm  thickness). 
Additionally,  the  thickness  of  the  midlayer  was  held  at  0.5mm  and  the  tri¬ 
material  configuration  was  given  a  temperature  increase  of  100°C.  The 
material  properties  for  the  six  cases  were  as  follows:  EA=EB=EC=  100  GPa, 
a^lOOxlOVC,  aB=  300xl0-VC,  and  Oc=  lOxlO’V’C. 

A.  NORMAL  STRESS  DISTRIBUTION 

Figure  5.2  contains  plots  of  the  normal  stresses  across  the  upper  and 
lower  interfaces  of  material  B  for  the  six  cases.  To  allow  comparison  of  the 
different  cases  the  abscissa  was  constructed  to  be  a  non-dimensional  location, 
with  the  distance  along  the  interface  being  divided  by  the  total  midlayer  length 
for  that  case. 

In  all  six  cases,  the  normal  stresses  remained  bearing  across  the  entire 
span  for  both  the  upper  and  lower  interfaces.  Cases  12,  16,  17,  and  18  (with 
midlayer  lengths  of  5mm,  9mm,  7.5mm,  and  2.5mm,  respectfully)  displayed 
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Case  17  LB=  7.5mm 


Case  12  LB=  5mm 


Case  18  LB=  2.5mm 


Case  19  LB=  1mm 


Case  20  LB=  0.5mm 

Figure  5.1  Tri-material  Configurations  For  Cases  12, 16, 17, 18, 19,  And 20 
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Figure  5.2  Normal  Stress  At  Interfeces  As  A  Function  Of  Non-dimensional 
Location  For  Cases  12, 16, 17, 18, 19,  And 20 
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similar  characteristics.  The  maximum  bearing  stresses  occurred  at  the  left 
and  right  ends  of  the  midlayer  with  the  minimum  being  developed  at  midspan. 
For  midlayer  lengths  of  9mm  and  7.5mm,  the  stress  along  the  midspan  was 
negligible.  As  the  length  was  further  decreased  to  5mm  and  2.5mm,  the 
midspan  normal  stresses  increased.  The  stresses  at  the  left  and  right  ends 
similarly  increased  but  to  a  much  lesser  extent  as  the  length  of  the  above 
mentioned  four  cases  decreased. 

Cases  19  and  20,  with  the  short  midlayer  lengths  of  1mm  and  0.5mm, 
displayed  distinctly  different  normal  stress  distributions  than  the  four  other 
cases.  The  left  ends  of  the  interfaces  had  relatively  low  bearing  stresses. 
Across  the  remaining  90%  of  the  span,  the  stress  distribution  was  quite 
uniform.  Figures  5.3  and  5.4  display  the  normal  stresses  acting  along  the 
upper  and  lower  interfaces  of  the  midlayer. 

B.  SHEAR  STRESS  DISTRIBUTION 

Plots  of  the  shear  stresses  along  the  upper  and  lower  interfaces  of  the 
midlayer  for  cases  12,  16,  17,  18,  19,  and  20  are  displayed  in  Figure  5.5.  Like 
the  normal  stress  distribution  plots,  the  abscissa  was  constructed  to  be  a  non- 
dimensional  location,  with  the  distance  along  the  interface  being  divided  by  the 
total  midlayer  length  for  that  case. 

All  six  cases  showed  similar  attributes  with  both  the  upper  and  lower 
shear  stresses  tending  to  axially  compress  the  midlayer.  Additionally,  the 
peak  stresses  occurred  at  the  extreme  left  and  right  ends  of  the  interfaces.  For 
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Material  B  Length  (mm)  I  5  I  9  I  7.5  I  2.5  I  1  1  0.5 

aA=  lOOxlO^C  aB=300xl0'6/°C  ac=  lOxlO^C  EA  =  EB=Ec=100GPa 
Figure  5.3  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  12, 16,  And  17 


Figure  5.5  Shear  Stress  At  Interlaces  As  A  Function  Of  Non-dimensional 
Location  For  Cases  12, 16, 17, 18, 19,  And  20 
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the  greatest  length  of  material  B,  case  16  with  1^=  9mm,  the  only  non- 
negligible  shear  stresses  appeared  at  the  ends  of  the  interfaces.  As  the 
midlayer  length  decreased,  the  extent  of  the  span  that  experienced  appreciable 
shear  stress  increased.  This  is  clearly  displayed  in  Figures  5.3  and  5.4. 

C.  SUMMARY 

The  effect  of  changing  the  midlayer  length  on  the  tri-material  configura¬ 
tion  was  shown  in  numerous  ways.  As  the  material  B  length  is  shortened,  the 
normal  stresses  tend  to  increase  and  become  uniformly  distributed  across  the 
span.  The  shear  stress  became  distributed  across  a  greater  extent  of  the 
interface  as  the  length  of  the  midlayer  lessened.  The  nature  of  the  normal 
stresses  was  always  bearing  and  were  much  less  than  the  shear  stresses 
generated.  As  the  length  of  the  midlayer  was  decreased  to  1mm,  shear  and 
normal  stresses  increased.  However,  further  reduction  in  length  resulted  in 
a  slight  decrease  of  both  stresses.  Except  for  the  shortest  length  where  it  was 
at  the  right  hand  edge,  the  maximum  shear  stresses  occurred  at  the  left  hand 
edge  for  the  upper  and  lower  interfaces.  The  maximum  normal  stresses  on 
both  interfaces  stayed  at  the  right  edge.  Finally,  the  maximum  upper  normal 
stresses  were  10  to  30  percent  greater  than  the  maximum  lower  normal 
stresses. 
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VI.  EFFECT  OF  NONUNIFORM  TEMPERATURE  DISTRIBUTION 


Chapters  III,  IV,  and  V  c<  <nsidered  the  tri-material  configuration  as  having 
a  uniform  temperature  increase  of  100°C  (AT  =  100°C).  In  this  chapter  the 
temperature  increases  of  the  three  materials  were  allowed  to  vary,  and  the 
results  were  studied  and  compared.  The  physical  dimensions  of  the  tri¬ 
material  configuration  used  in  this  chapter  are  shown  in  Figure  6. 1.  Materials 
A  and  C  were  both  given  a  length  of  10mm  and  a  height  of  1mm.  Material  B 
was  given  a  length  of  2.5mm  and  a  height  of  0.1mm.  The  material  properties 
supplied  for  all  cases  studied  were  as  follows:  EA  =  EB  =  Ec  =  100  GPa,  aA  = 
lOOxlOVC,  <xB  =  300xl0VC,  and  ac  =  lOxlOVC. 


Figure  6.1  Tri-material  Configuration  For  Cases  21, 22, 23,  And  24 


Four  different  cases  were  examined.  In  case  2 1,  all  three  materials  were 
allowed  a  temperature  rise  of  100°C  (ATa  =  ATB  =  ATC  =  100°C).  A  tempera¬ 
ture  increase  of  200°C  for  materials  A  B,  and  C  (ATa  =  ATB  =  ATC  =  200°C) 
was  studied  in  case  22.  For  case  23,  material  A  was  assigned  a  temperature 
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rise  of  200°C  while  materials  B  and  C  were  permitted  100°C  temperature 
increases  (ATa  =  200°C,  ATB  =  ATC  =  100°C).  FinaLy,  for  case  24,  material  A 
was  allowed  to  have  a  linearly  decreasing  temperature  profile.  The  top  layer 
of  material  A  was  set  at  a  temperature  increase  of  200°C  and  the  temperature 
rise  was  then  lessened  until  at  the  interface  of  materials  A  and  B  the 
temperature  increase  was  100°C.  This  profile  is  shown  in  Figure  6.2. 
Materials  B  and  C  were  allowed  100°C  rises  in  temperature  (ATa  =  linear,  ATB 
=  ATC  =  100°C). 

200  °C 

150  °C 

125  °C 
113  °C 
106  °C 
100  °C 

Figure  6.2  Material  A  Linear  Temperature  Distribution  For  Case  24 
A.  NORMAL  STRESS  DISTRIBUTION 

The  normal  stresses  acting  along  the  upper  and  lower  interfaces  of  the 
midlayer  for  the  four  cases  are  shown  in  Figure  6.3.  Similar  characteristics 
were  displayed  for  both  the  upper  and  lower  normal  stress  distributions.  The 
maximum  stresses  occurred  at  the  far  left  and  right  ends  of  the  interfaces  and 
were  always  bearing.  For  the  majority  of  the  midspan  the  normal  stresses 
were  minimal. 
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Figure  6.3  Normal  Stresses  At  The  Interfeces  For  Cases  21, 22, 23,  And  24 
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A  comparison  between  case  21  and  case  22  revealed  an  interesting  point. 
Case  22,  with  a  tri-material  configuration  temperature  rise  of  200°C,  resulted 
in  normal  stresses  which  were  twice  the  value  of  the  normal  stresses  of  case 
21  (ATa  =  ATb  =  ATC  =  100°C).  This  showed  the  direct  relationship  between 
system  temperature  increase  and  the  resulting  increase  in  normal  stresses. 

When  only  material  A  was  given  a  larger  temperature  rise,  as  in  case  23 
(ATa  =  200°C,  ATB  =  ATC  =  100°C),  the  change  in  system  behavior  for  normal 
stresses  when  compared  to  a  uniform  temperature  field  (case  2 1)  was  minimal. 
Along  the  upper  interface  the  response  was  nearly  identical,  with  the  only 
difference  being  the  marginally  larger  bearing  stress  at  the  extreme  right  end 
for  case  23.  Across  the  lower  interface,  case  23  normal  stress  values  were 
continuously  slightly  less  than  those  of  case  21. 

When  the  temperature  rise  in  material  A  was  allowed  to  linearly  decrease 
through  its  thickness,  as  in  case  24,  the  resulting  normal  stress  distributions 
were  the  lowest  of  the  four  cases  studied.  The  maximum  stresses  developed 
were  approximately  one-half  of  the  maximum  for  the  uniform  temperature  field 
case. 

Figures  6.4  and  6.5  display  the  normal  stress  distributions  acting  on  the 
upper  and  lower  interfaces  of  the  midlayer  for  the  four  cases. 

B.  SHEAR  STRESS  DISTRIBUTION 

Figure  6.6  contains  the  plots  of  the  shear  stress  distributions  across  the 
upper  and  lower  interfaces  for  cases  21,  22,  23,  and  24.  The  resulting  curves 


60 


336  MPa 


« 


3  1 

+ 

1 

n 

H 

♦  ♦ 

CM 

H 

c3 

§ 

8 

8 

< 

-S 

^-4 

O 

h3 

m 

CM 

C 

1 

8 

8 

O 

CM 

CM 

w 

cn 

< 

200 

i 

8 

CM 

u 

CM 

w 

Xfl 

8 

8 

8 

< 

H 

o 

^■s 

/“N 

o 

O 

O 

< 

CQ 

c 

H 

H 

H 

< 

<3 

< 

61 


Figure  6.4  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  21  And  22 
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Figure  6.5  Shear  And  Normal  Stresses  On  The  Midlayer  For  C  ases  23  And  24 


LOWER  SHEAR  STRESS  (MPa)  UPPER  SHEAR  STRESS  (MPa) 


Figure  6.6  Shear  Stresses  At  The  Interfaces  For  Cases  21, 22, 23,  And  24 
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were  similar  for  all  four  cases.  The  maximum  shear  stresses  occurred  at  the 
extreme  right  and  left  ends  of  the  interfaces.  Along  the  majority  of  the  length 
of  span  the  shear  stresses  fluctuated  around  zero.  This  fluctuation  may  have 
been  due  to  a  breakdown  in  the  finite  element  program  in  dealing  with  the 
very  small  numerical  values  present. 

As  with  the  normal  stresses,  a  comparison  of  case  21  (ATa  =  ATB=  ATC  = 
100°C)  with  case  23  (ATa  =  ATB  =  ATC  =  200°C)  shows  that  the  shear  stresses 
at  every  point  for  the  higher  temperature  increase  were  twice  the  value  of 
those  stresses  for  the  lower  temperature  rise.  This  also  displayed  the  direct 
relationship  between  system  temperature  increase  and  the  resulting  increase 
in  shear  stresses. 

Case  23  (ATa  =  200°C,  ATB  =  ATC  =  100°C)  and  case  24  (ATa  =  Unear,  ATB 
=  ATC  =  100°C)  produced  shear  stress  distribution  values  very  close  to  case  2 1. 
This  suggests  that  for  the  particular  tri-material  configuration  chosen  for  the 
chapter,  increasing  the  temperature  rise  for  one  material  does  not  appreciably 
change  the  system's  shear  stress  response. 

Displays  of  the  shear  stress  distributions  acting  on  the  upper  and  lower 
interfaces  for  the  four  cases  are  shown  in  Figures  6.4  and  6.5. 

C.  SUMMARY 

A  study  of  the  effects  of  a  nonuniform  temperature  field  produced  a 
number  of  interesting  results.  A  doubling  of  the  system  temperature  rise 
resulted  in  a  doubUng  of  the  normal  and  shear  stresses.  Normal  stresses 
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remained  bearing  in  different  uniform  and  nonuniform  temperature  fields. 
The  normal  and  shear  stress  distributions  maintained  the  same  general  shapes 
regardless  of  temperature  distribution.  Increasing  the  temperature  rise  in  just 
one  material  resulted  in  minimal  change  in  system  response.  Finally,  a 
temperature  rise  which  was  allowed  to  linearly  decrease  through  the  thickness 
of  one  material  resulted  in  the  lowest  normal  stress  distributions. 
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VII.  SUMMARY  AND  CONCLUSIONS 


The  results  of  this  investigation  into  the  behavior  of  a  tri-material  medium 
subjected  to  a  temperature  field  leads  to  several  conclusions.  These 
conclusions  are  broken  down  into  the  following  topics:  effect  of  material 
properties,  effect  of  thickness  of  midlayer,  effect  of  length  of  midlayer,  and 
effect  of  nonuniform  temperature  field. 

A.  EFFECT  OF  MATERIAL  PROPERTIES 

When  both  the  temperature  rise  and  the  coefficients  of  thermal  expansion 
were  kept  uniform  throughout  the  system  and  only  Young's  Modulus  was 
allowed  to  vary,  no  normal  or  shear  stresses  were  developed. 

Changing  the  coefficient  of  thermal  expansion  produced  several  results. 
As  the  mismatch  between  the  coefficients  of  thermal  expansion  of  the  materials 
became  greater,  the  normal  and  shear  stress  magnitudes  along  the  upper  and 
lower  interfaces  became  greater.  The  shear  stresses  developed  along  the 
interfaces  were  much  greater  than  the  normal  stresses.  Depending  on  the 
values  of  the  coefficients  of  thermal  expansion,  the  normal  stresses  were  either 
bearing  or  peeling.  Finally,  along  the  upper  and  lower  interfaces  the  normal 
and  shear  stresses  were  very  local  in  nature.  Only  in  the  immediate  area  of 
the  left  and  right  corners  of  the  interface  were  the  stresses  significant. 
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B.  EFFECT  OF  MIDLAYER  THICKNESS 


For  the  particular  tri-material  configuration  studied  the  normal  stresses 
remained  bearing  as  the  midlayer  became  thinner.  The  maximum  normal  and 
shear  stresses  decreased  with  decreasing  material  B  thickness.  As  the 
midlayer  was  made  thinner,  the  normal  and  shear  stress  distributions  became 
more  and  more  local  in  nature.  It  is  believed  that  as  the  midlayer  becomes 
extremely  thin,  the  only  stresses  of  any  consequence  would  exist  at  the 
extreme  left  and  right  corners  of  the  interfaces.  The  shear  stresses  generated 
were  two  to  three  times  greater  than  the  normal  stresses.  As  the  thickness  of 
material  B  was  decreased,  the  location  of  the  maximum  shear  stress  along  the 
lower  interface  moved  from  the  inner  edge  to  the  outer  edge.  However,  along 
the  upper  interface  the  maximum  shear  stress  remained  fixed  at  the  inner 
edge.  Lastly,  the  maximum  normal  stress  stayed  at  the  left  edge  for  both  the 
upper  and  lower  interfaces. 

C.  EFFECT  OF  MIDLAYER  LENGTH 

With  the  tri-material  configuration  used  the  resulting  shear  stresses  were 
much  greater  than  the  normal  stresses.  The  normal  stresses  developed  were 
always  bearing.  As  midlayer  length  was  decreased  to  1mm,  the  normal  and 
shear  stresses  increased.  Further  shortening  of  the  length  of  the  midlayer 
resulted  in  slight  decrease  of  the  shear  and  normal  stresses.  On  both 
interfaces  the  maximum  shear  stress  occurred  at  the  left  hand  edge.  The  only 
exception  was  the  shortest  length,  where  the  maximum  shear  stress  appeared 
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at  the  right  hand  edge.  The  maximum  normal  stress  was  experienced  at  the 
right  hand  edge  for  both  the  upper  and  lower  interfaces.  Additionally,  the 
upper  maximum  normal  stress  was  approximately  10  to  30  percent  larger  than 
the  lower  maximum  normal  stress. 

D.  EFFECT  OF  NONUNIFORM  TEMPERATURE  DISTRIBUTION 

When  the  temperature  rise  for  all  three  materials  of  the  system  was 
doubled,  the  result  was  a  doubling  of  the  normal  and  shear  stresses.  The 
normal  stresses  remained  bearing  and  the  shear  and  normal  stress  distribu¬ 
tions  maintained  the  same  general  shapes  regardless  of  the  temperature 
distribution.  Increasing  the  temperature  rise  in  just  the  upper  layer  resulted 
in  minimal  change  in  system  response.  Finally,  a  linearly  decreasing 
temperature  rise  gave  the  lowest  magnitude  normal  stress  distributions. 

E.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

A  number  of  different  areas  for  further  research  in  this  topic  exist. 
Validate  the  tri-material  response  to  nonuniform  temperature  fields  by  using 
a  larger  number  of  system  configurations.  Alter  the  finite  element  code  to 
allow  study  of  the  system  response  to  transient  temperature  fields.  Refine  the 
tri-material  configuration  to  more  closely  match  electronic  component 
packaging.  Finally,  investigate  the  effect  of  thermal  cycling  on  fatigue. 
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APPENDIX  A 
FEM  PROGRAM 


FEM  PROGRAM  FOR  WELD  PROBLEM 


COMPILE  AND  RUN  THIS  CODE  WITH  THE  SINGLE  EXEC  COMMAND:  WELD 


NEL  THE  NUMBER  OF  ELEMENTS 

NODES  THE  NUMBER  OF  SYSTEM  NODES 

NDOF  THE  NUMBER  OF  SYSTEM  DEGREES  OF  FREEDOM 

NBC  THE  NUMBER  OF  ESSENTIAL  BOUNDARY  CONDITIONS 

IBC  THE  SUPPRESSED  DEGREES  OF  FREEDOM 

NFC  THE  NUMBER  OF  APPLIED  FORCES 

IFC  THE  SYSTEM  NODES  WITH  APPLIED  LOADS 

I CORK  THE  CORRESPONDENCE  BETWEEN  LOCAL  AND  SYSTEM  DOF 

NMAT  NUMBER  OF  DIFFERENT  MATERIALS 

IMAT  A  PARTICULAR  MATERIAL 

ITYPE  -  0  FOR  B.C. 'S  AT  LEFT  END  ONLY, 

-  1  FOR  B.C. 'S  AT  BOTH  LEFT  AND  RIGHT  ENDS 
E  YOUNGS  MODULUS  OF  IMAT  MATERIAL 

G  SHEAR  MODULUS  OF  IMAT  MATERIAL 

X  THE  X  COORDINATE  OF  A  NODE 

HITE  THE  HEIGHT  OF  AN  ELEMENT 

ZLEN  THE  LENGTH  OF  AN  ELEMENT 

EK  THE  6  BY  6  ELEMENT  STIFNESS  MATRIX 

SYSMAT  THE  NDOF  BY  NDOF  SYSTEM  STIFFNESS  MATRIX 
FLO  AD  THE  MAGNITUDE  OF  AN  APPLIED  LOAD 

FORCE  THE  SYSTEM  FORCE  VECTOR 

U  THE  SYSTEM  DISPLACEMENT  VECTOR 

DELTMP  THE  TEMPERATURE  CHANGE 

ALPHA  THERMAL  COEFFICIENT  OF  EXPANSION 
NNSCOL  THE  NUMBER  OF  NONSOLDERED  COLUMNS 
NSCOL  THE  NUMBER  OF  SOLDERED  COLUMNS 
N1ROW  THE  NUMBER  OF  ROWS  IN  THE  CHIP 

N3ROW  THE  NUMBER  OF  ROWS  IN  THE  SUBSTRATE 

PLSTRS  THE  PEELING  STRESS  AT  THE  NODE  TO  THE  RIGHT  OF  THE 
COLUMN  NUMBER 

N3STRT  ELEMENT  NUMBER  BEFORE  START  OF  LOWER  SOLDER 

INCLUDE  ' WELDCOMN ' 

CHARACTER* 20  NAME 
COMMON  / WORKSP/ RWXSP 
REAL  RWKSP( 137047) 

CALL  IWKIN(137047) 

READ  *,  NAME 

NOTE;  INPUT  THE  NAME  OF  THE  FILE  TO  BE  OPENED  AS  FOLLOWS: 

'/FN  FT' 

THE  FILENAME  MUST  BE  PRECEDED  WITH  A  '/  AND  THE  FILETYPE  MUST  BE 

FOLLOWED  BY  A  THE  ONLY  SPACE  IS  BETWEEN  FN  AND  FT. 

OPEN (9,  FILE-  NAME  ) 

INPUT  THE  FEM  PARAMETERS 

READ(9,*)  NEL,  NODES,  NMAT,  NDOF,  NBC,  NFC,  ITYPE 
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*  READ(9,*)  (IBC(I),  X  -  1,NBC) 

DO  50  XB  -  1,NBC 
50  CONTINUE 

Xr  (NFC  .NE.  0)  THEN 

READ {  9 ,  * )  (IFC(X) ,  FLOAD(I),  X  -  1,NFC) 

DO  70  IF  •  1,NFC 
70  CONTINUE 
ELSE 

ENDXF 

DO  80  XH  -  l.NMAT 

READ (9,81)  E(XM) ,  POISS(IM),  ALPHA ( IM) ,  HXTE(XM) 

80  CONTINUE 

81  FORMAT  (4F10.7) 


READ (9, 22)  (X(IN),  IN  -  1, NODES) 

DO  999  ILL  -  1,  NEL 

READ (9,20)  XL, IMAT(ILL) , (ICORR(ILL, IN) , IN  -  1,6),  INODE (ILL, 1) , 

&  INODE (ILL, 2 ) 

999  CONTINUE 
20  FORMAT (1014) 

22  FORMAT(7F10 . 5) 

READ (9, *)  DELTMP 

READ (9, *)  (TEMP (I) ,1  -  l.NMAT+l) 

WRITE (14 , *)  '  1 

WRITE (14,*)  'DELTMP  -  ',  DELTMP 

* 

*  CALCULATE  THE  SHEAR  MODULI I 

• 

WRITE (14 , *)  '  ' 

WRITE ( 14 , *)  ' I MAT/ ,  E,  6  MODULI!,  ALFA,  POISSON  AND  HEIGHT  FOLLOW' 
WRITE (14 ,  *)  '  ' 

DO  100  IM  -  1,  NMAT 

G  (IM)  -  5. *E(IM) / (12 . * (l+POISS(IM) ) ) 

WRITE ( 14 , 105)  IM,  E(IM) ,  G(IM) ,  ALPHA (IM) ,  POISS(IM),  HITE(IM) 

100  CONTINUE 

105  FORMAT(3X,  13 , 2X, 5E12 . 5) 

* 

WRITE (14 , *)  '  ' 


DO  350  I  -  l,NDOF 
FORCE ( I )  -  0.0 
DO  340  J  -  1 , NDOF 
SYSMAT(I, J)  -0.0 
340  CONTINUE 
350  CONTINUE 

* 

•  THE  ELEMENT  DO  LOOP  TO  500  FORMS  THE  SYSTEM  MATRIX  AND  VECTOR 

* 

DO  500  I EL  -  1 , NEL 
IM  -  IMAT(IEL) 

* 

*  FORM  THE  ELEMENT  STIFFNESS  MATRICES 

* 

IL  -  INODE (IEL, 1) 

XR  -  INODE (IEL, 2) 
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p 


* 


* 


* 


. ZLEN(XEL)  -  X(IR)  -  X(XL) 

H  -  HITE(XM) 

EE  -  E(XM) 

CG  -  G(XM) 

A  -  GG*ZLEN ( IEL) / ( 4 • *H) 

B  -  GG*H/ZLEN(IEL) 

C  -  EE*H/ ( 6 . *ZLEN ( IEL)  ) 

D  -  0 . 5*GG 
AP2C  -  A  ♦  2 . *C 
ZMAPC  -  -  A  ♦  C 
AM2C  -  A  -  2.*C 
ZMAMC  -  -  A  -  C 

EK(1, 1)  -  AP2C 
EK ( 1 , 2 )  -  ZMAPC 
EK  ( 1 , 3 )  -  D 
EK ( 1 , 4 )  -  AM2C 
EK ( 1 , 5)  -  ZMAMC 
EK ( 1 , 6 )  -  -D 
EK (2 , 1)  -  ZMAPC 
EK ( 2 , 2 )  -  AP2C 
EK (2,3)  -  -D 
EK  (2 , 4 )  -  ZMAMC 
EK (2,5)  -  AM2C 
EK (2 , 6)  -  0 
EK  ( 3 , 1 )  -  0 
EK ( 3 , 2 )  -  -D 
EK (3,3)  -  B 
EK (3,4)  -  D 
EK (3 , 5)  -  -D 
EK ( 3 , 6)  -  -B 
EK  (4 , 1)  -  AM2C 
EK (4 , 2 )  -  ZMAMC 
EK(4 , 3 )  -  D 
EK(4,4)  -  AP2C 
EK (4,5)  -  ZMAPC 
EK (4,6)  -  -D 
EK (5,1)  -  ZMAMC 
EK(5, 2)  -  AM2t 
EK ( 5 , 3 )  -  -D 
EK ( 5 , 4 )  -  ZMAPC 
EK(5,5)  -  AP2C 
EK (5,6)  -  0 
EK  ( 6 , 1)  -  -D 
EK(6,2)  -  D 
EK ( 6 , 3 )  -  -B 
EK (6 , 4)  -  -D 
EK ( 6 , 5 )  -  0 
EK  (6 , 6)  -  B 

ASSEMBLE  ELEMENT  MATRICES  INTO  GLOBAL  STXFFNES  MATRIX 


DO  480  ID  -  1,6 
IN  -  ICORR(XEL, ID) 

DO  480  JD  -  1,6 
JN  -  ICORR(IEL.JD) 

SYSMAT(XN.JN)  -  SYSMAT(IN,JN)  +  EK(ID,JD) 
480  CONTINUE 
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<  500  CONTINUE 

« 

*  CONSTRUCT  THE  FORCE  VECTOR  DUE  TO  TEMPERATURE  CHANGE 

* 

WRITE (14,*)  •  • 

IF  (DELTMP  .NE.  0.0)  THEN 
DO  550  I EL  -  1,  NEL 
IM  -  IMAT(IEL) 

CHI2  -  E ( IM) ‘ALPHA ( IM) * (2 . *TEMP (IM) +TEMP ( IM+1) )*HIT£(IM)/6. 
CHI1  -  E(IM)  *ALPHA(IM)  *  (TEMP(IM)  +2  .  *TEMP(IM+1) )  *HITE(IM)  /  6. 

11-  ICORR(IEL, 1) 

12-  ICORR(IEL, 2) 

14-  I CORK ( IEL, 4 ) 

15-  ICORR(IEL,5) 

FORCE (II)  -  FORCE (II)  -  CHI1 
FORCE (12)  -  FORCE (12)  -  CHI2 
FORCE (14)  -  FORCE (14)  +  CHI1 
FORCE (IS)  -  FORCE (15)  +  CHI2 
550  CONTINUE 
ELSE 
Et'DIF 

* 

*  MODIFY  STIFFNESS  MATRIX  AND  FORCE  VECTOR  FOR  BOUNDARY  CONDITIONS 

* 

DO  610  1  -  l.NBC 
II  -  IBC(I) 

DO  600  J  -  l.NDOF 
SYSMAT (II , J)  -  0.0 
600  CONTINUE 

SYSMAT (II, II)  -  1. 

FORCE (II)  -  0.0 
610  CONTINUE 

* 

IF ( ITYPE  -.EQ.  1)  THEN 
READ (9,*)  ILASTTOP 
WRITE (14 , *)  '  ' 

WRITE (14 , *)  'B.C.  S  AT  BOTH  ENDS;  ILASTTOP  -  ILASTTOP 
DO  620  IN  -  l,NDOF 
SYSMAT (ILASTTOP, IN)  -  0.0 
620  CONTINUE 

SYSMAT (ILASTTOP, ILASTTOP)  -  1. 

SYSMAT (ILASTTOP , NDOF)  -  -1. 

FORCE (ILASTTOP)  -  0.0 
ELSE 
ENDIF 

* 

IF  (NFC  .NE.  0)  THEN 
DO  630  IF  -  1,NFC 
II  -  IFC(IF) 

FORCE (II)  -  FORCE (II)  +  FLOAD(IF) 

630  CONTINUE 
ELSE 
ENDIF 

* 

*  SOLVE  THE  SYSTEM  OF  EQUATIONS  FOR  THE  DISPLACEMENTS 

* 

CALL  DLSARG (NDOF,  SYSMAT,  NPARAM,  FORCE,  1,  U) 

* 

*  PRINT  THE  FORCE  AND  DISPLACEMENT  VECTORS 
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WRITE (14,*)  «  ' 

«  WRITE (14,*)  '  rORCE  AND  DISPLACEMENT  VECTORS  FOLLOW' 

WRITE(14,*)  •  ' 

DO  680  I DOF  -  1 , NDOF 

WRITE ( 14 , 641)  IDOF ,  FORCE (I DOF ) ,  U(IDOF) 

680  CONTINUE 

641  FORMAT(2X,  14,  4X,  E12.5,  10X,  E12.5) 

* 

*  CALCULATE  THE  STRESSES 

* 

DO  700  I EL  «  1,  NEL 

11  -  ICORR(IEL, 1) 

12  •  ICORR(I£L,2) 

13  -  ICORR(IEL,3) 

14  «  I CORK (I EL, 4) 

15  -  ICORR(IEL, S) 

16  -  ICORR(IEL, 6) 

U1BOT  -  U(I1) 

U1TOP  -  U(I2) 

VI  -  U(I3) 

U2BOT  -  U( 14 ) 

U2TOP  ■  U(I5) 

V2  -  U(I6) 

IM  -  IMAT(IEL) 

EE  -  E(IM) 

GG  -  G(IM) 

SIGMABOT (IEL) “  EE* ( (U2BOT  -  U1B0T) /ZLEN(IEL)  -  ALPHA (IM)* 

6  TEMP(IM+1) ) 

SIGMA TOP (IEL) -  EE*((U2TOP  -  U1TOP) /ZLEN(IEL)  -  ALPHA ( IM) *TEMP(IM) ) 
TAUELEM(IEL)  -  GG*(  ( (U1TOP  -U1BOT)  +■  (U2TOP  -  U2BOT)  )  / 

6  (2 . *HITE(IM) )  +  (V2  -  VI) /ZLEN(IEL)  ) 

700  CONTINUE 

* 

WRITE(14 ,*)  •  • 

WRITE (14, *)  '  STRESSES  FOLLOW' 

WRITE (14, *)  '  ' 

WRITE (14,*)  '  SGMABOT  SGMATOP  TAUELEM  ' 

* 

PSI  -  E ( 1 ) *ALPHA ( 1 ) *DELTMP 
DO  800  IEL  -  1,NEL 
IM  -  IMAT (IEL) 

SIGBTNDL  -  SIGMABOT (IEL) /PSI 
SIGTPNDL  -  SIGMATOP(IEL) /PSI 
TAUNDL  -  TAUELEM (I EL) /PSI 

WRITE (14, 801)  IEL,  SIGMABOT ( IEL) ,  SIGMATOP(IEL) ,  TAUELEM ( IEL) 

800  CONTINUE 

801  FORMAT ( (2X,I4,6(4X,G12.S))) 

* 

WRITE ( 14 , *)  '  ' 

* 

READ(9, *}  NNSCOL,  NSCOL,  N1ROW,  N3ROW 

* 

NCOL  »  NNSCOL  +  NSCOL 

* 

DO  1300  IEL  -  6,  70 
IM  -  IMAT (I EL) 

PLFCE(IEL)  -  TAUELEM (I EL) *HITE(IM) 

1300  CONTINUE 
II  -  6 


73 


*  o 


DO  1500  I COL  -  1,8 
*  I ELF  -  6  ♦  (ICOL  -  1) 

I ELL  »  XELF  +  64 
PLF  -  0.0 

DO  1400  XEL  -  XELF, I ELL, 14 

PLF  -  PLF  +  PLFCE(IEL+1)  -  PLFCZ(IEL) 

1400  CONTINUE 

UPLSTR(XCOL)  -  2 . *PLF/ (X(II+2)  -  X(II)) 

WRITE (14,31)  ICOL,  UPLSTR(ICOL) 

XI  -  XI  ♦  1 
1500  CONTINUE 

31  FORMAT (IX, 12, '  UPPER  PEEL  STRESS  -  ',G12.5) 

• 

WRITE (14,*)  '  ' 

DO  1600  I EL  -  121,  197 
IM  -  IMAT(IEL) 

PLFCE(IEL)  -  TAUELEM(IEL) *HITE(IM) 

1600  CONTINUE 
II  -  6 

DO  1800  ICOL  -  1,8 
I ELF  -  121  ♦  (ICOL  -  1) 

I ELL  -  XELF  +  76 
PLF  -  0.0 

DO  1700  XEL  -  IELF , I ELL, 17 

PLF  -  PLF  +  PLFCE (IEL+1)  -  PLFCE(IEL) 

1700  CONTINUE 

BPLSTR(ICOL)  «(2.*PLF/(X(II+2)  -  X (II) ) ) * (-1 . 0) 

WRITE (14,32)  ICOL, BPLSTR (ICOL) 

II  -  II  +  1 
1800  CONTINUE 

32  FORMAT (IX, 12, '  LOWER  PEEL  STRESS  -  \G12.5) 

* 

*  COMPARE  MOMENTS  DUE  TO  BENDING  AND  PEEL  STRESSES 

* 

AVEFCE(l)  -  0. 5* (SIGMATOP(l)  +  SIGMABOT(l) ) *HITE(1) 
ZMOMARM ( 1 )  -  .0562  -  .014 

AVEFCE ( 2 )  -  0.5* (SIGMATOP ( 15)  +  SIGMABOT(15) ) *HITE(2) 
ZMOMARM  ( 2 )  -  .0562  -  .028  -  .007 

AVEFCE ( 3 )  -  0.5* (SIGMATOP (29)  +  SIGMABOT(29) ) *HITE(3) 
ZMOMARM ( 3 )  -  .0562  -  .028  -  .014  -  .0035 

AVEFCE (4)  -  0.5* (SIGMATOP (4 3)  +  SIGMABOT(43) ) *HITE(4) 
ZMOMARM ( 4 )  -  .0035  *  .00175 

AVEFCE (5)  -  0.5* (SIGMATOP (57)  ♦  SIGMABOT(57) ) *HITE(4) 
ZMOMARM (5)  -  .00175 
ZMOMBDG  «  0.0 
DO  1900  1-1,5 

ZMOMBDG  -  ZMOMBDG  +  ZMOMARM (I) * AVEFCE (I) 

1900  CONTINUE 

• 

ZMOMPEEL  -  0.0 
DO  2000  I  -  1,8 
II  -  7  ♦  I 
JJ  -  II  -  2 
KX  »  II  -  1 

W  -  0 . 5* (X (II)  -  X(JJ)) 

Y  -  0.25*(X(JJ)  +  X(II)  +2 . *X (KX) ) 

ZMOMPEEL  -  ZMOMPEEL  +  W*Y*UPLSTR(I) 

2000  CONTINUE 

WRITE (14,*)  'MOMPEEL  -  ',  ZMOMPEEL 
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<  2M0MBDG  -  0.0 
IM  -  0 
Y2  -  0.056 
DO  3000  I  -  1,57,14 
IH  -  XMAT(I) 

Y1  -  Y2  -  HITE(IM) 

AOV3  -  (SIGMATOP(X) -SIGMABOT(I) ) / (3 . *HITE(IM) ) 

BOV2  -  0.5*(SICHATOP(I)-(Y2/HITE(IM) ) * (SIGMATOP(I) -SIGMABOT(I) ) ) 
2HOMBDG  -  2HOMBDG  +  AOV3* ( Y2**3-Y1**3) +BOV2* ( Y2**2-Y1**2 ) 

Y2  -  Y1 
3000  CONTINUE 

SXGMAT  ■  1.052 
SXGMAB  -  2.3012 
Y2  •  .001 
Y1  -  .0005 

B  -  SIGMAT  -  (Y2/ . 0005) * (SXGMAT-SIGKAB) 

A  -  (SXGMAT-SXGMAB) / ( . 0005) 

STOP 

END 
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APPENDIX  B 
PROGRAM  EXEC 

FORTVS2  WELD  (AD(DBL4)  ICA) 

FILEDEF  14  DISK  WELD  OUTPUT  (RECFM  FB  LRECL  102  BLOCK  10200) 
LOAD  WELD  (START 
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APPENDIX  C 
INPUT  FILE 


i 


fOO  18  15 
1  2  3  4  5 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09 
100.E+09. 
0.0 
5.50 
10. 


261  14  0  0 
6  7  146  164  165  166 
P.3  100.E-06 


.3 

.3 

.3 

.3 

.3 

.3 

2<  3 

.3 
.3 
.3 
.3 
.3 
.3 
..3 

1.33333 

6.0 


100.E-06 

100.E-06 

100.E-06 

100.E-06 

300.E-06 

300.E-06 

300.E-06 

300.E-06 

300.E-06 

10.E-06 

10.E-06 

10.E-06 

10.E-06 

10.E-06 

2.66667 

7.0 


10.000001  10.000002 


167  168  169 
.500 
.250 
.125 
.0625 
.0625 
.010 
.  020 
>3 

.020 
.010 
.0625 
.0625 
.125 
.250 
•  500, _ 
4.0 
8.0 

10.000003 


1 

1 

3 

1 

2 

10 

8 

9 

1 

2 

2 

1 

10 

8 

9 

17 

IS 

16 

2 

3 

3 

1 

17 

15 

16 

24 

22 

23 

3 

4 

4 

1 

24 

22 

23 

31 

29 

30 

4 

5 

5 

1 

31 

29 

30 

38 

36 

37 

5 

6 

6 

1 

38 

36 

37 

45 

43 

44 

6 

7 

7 

1 

45 

43 

44 

52 

SO 

51 

7 

8 

8 

1 

52 

50 

51 

59 

57 

58 

8 

9 

9 

1 

59 

57 

58 

66 

64 

65 

9 

10 

10 

1 

66 

64 

65 

73 

71 

72 

10 

11 

11 

1 

73 

71 

72 

80 

78 

79 

11 

12 

12 

1 

80 

78 

79 

87 

85 

86 

12 

13 

13 

1 

87 

85 

86 

94 

92 

93 

13 

14 

14 

1 

94 

92 

93 

101 

99 

100 

14 

15 

15 

2 

4 

3 

2 

11 

10 

9 

1 

2 

16 

2 

11 

10 

9 

18 

17 

16 

2 

3 

17 

2 

18 

17 

16 

25 

24 

23 

3 

4 

18 

2 

25 

24 

23 

32 

31 

30 

4 

5 

19 

2 

32 

31 

30 

39 

38 

37 

5 

6 

20 

2 

39 

38 

37 

46 

45 

44 

6 

7 

21 

2 

46 

45 

44 

53 

52 

51 

7 

8 

22 

2 

53 

52 

51 

60 

59 

58 

8 

9 

23 

2 

60 

59 

58 

67 

66 

65 

9 

10 

24 

2 

67 

66 

65 

74 

73 

72 

10 

11 

25 

2 

74 

73 

72 

81 

80 

79 

11 

12 

26 

2 

81 

80 

79 

88 

87 

86 

12 

13 

27 

2 

88 

87 

86 

95 

94 

93 

13 

14 

28 

2 

95 

94 

93 

102 

101 

100 

14 

15 

29 

3 

5 

4 

2 

12 

11 

9 

1 

2 

30 

3 

12 

11 

9 

19 

18 

16 

2 

3 

31 

3 

19 

18 

16 

26 

25 

23 

3 

4 

32 

3 

26 

25 

23 

33 

32 

30 

4 

5 

33 

3 

33 

32 

30 

40 

39 

37 

5 

6 

34 

3 

40 

39 

37 

47 

46 

44 

6 

7 

35 

3 

47 

46 

44 

54 

53 

51 

7 

8 

36 

3 

54 

53 

51 

61 

60 

58 

8 

9 

37 

3 

61 

60 

58 

68 

67 

65 

9 

10 

38 

3 

68 

67 

65 

75 

74 

72 

10 

11 

39 

3 

75 

74 

72 

82 

81 

79 

11 

12 

40 

3 

82 

81 

79 

89 

88 

86 

12 

13 

4 


4.66667  5.0  5.25 

9.0  9.50  9.750 


1.  Young’s  Modulus 

2.  Poisson's  Ratio 

3.  Coefficient  of 
Thermal  Expansion 

4.  Row  Height 

5.  Column  Width 
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< 


41  3  89  88  86  96  95  93  13  14 

42  3  96  95  93  103  102  100  14  15 

43  4  6  5  2  13  12  9  1  2 

44  4  13  12  9  20  19  16  2  3 

45  4  20  19  16  27  26  23  3  4 

46  4  27  26  23  34  33  30  4  5 

47  4  34  33  30  41  40  37  5  6 

48  4  41  40  37  48  47  44  6  7 

49  4  48  47  44  55  54  51  7  8 

50  4  55  54  51  62  61  58  8  9 

51  4  62  61  58  69  68  65  9  10 

52  4  69  68  65  76  75  72  10  11 

53  4  76  75  72  83  82  79  11  12 

54  4  83  82  79  90  89  86  12  13 

55  4  90  89  86  97  96  93  13  14 

56  4  97  96  93  104  103  100  14  15 

57  5  7  6  2  14  13  9  1  2 

58  5  14  13  9  21  20  16  2  3 

59  5  21  20  16  28  27  23  3  4 

60  5  28  27  23  35  34  30  4  5 

61  5  35  34  30  42  41  37  5  6 

62  5  42  41  37  49  48  44  6  7 

63  5  49  48  44  56  55  51  7  8 

64  5  56  55  51  63  62  58  8  9 

65  5  63  62  58  70  69  65  9  10 

66  5  70  69  65  77  76  72  10  11 

67  5  77  76  72  84  83  79  11  12 

68  5  84  83  79  91  90  86  12  13 

69  5  91  90  86  98  97  93  13  14 

70  5  98  97  93  105  104  100  14  15 

71  6  106  42  37  107  49  44  6  7 

72  6  107  49  44  108  56  51  7  8 

73  6  108  56  51  109  63  58  8  9 

74  6  109  63  58  110  70  65  9  10 

75  6  110  70  65  111  77  72  10  11 

76  6  111  77  72  112  84  79  11  12 

77  6  112  84  79  113  91  86  12  13 

78  6  113  91  86  114  98  93  13  14 

79  6  114  98  93  115  105  100  14  15 

80  7  116  106  37  117  107  44  6  7 

81  7  117  107  44  118  108  51  7  8 

82  7  118  108  51  119  109  58  8  9 

83  7  119  109  58  120  110  65  9  10 

84  7  120  110  65  121  111  72  10  11 

85  7  121  111  72  122  112  79  11  12 

86  7  122  112  79  123  113  86  12  13 

87  7  123  113  86  124  114  93  13  14 

88  7  124  114  93  125  115  100  14  15 

89  8  126  116  37  127  117  44  6  7 

90  8  127  117  44  128  118  51  7  8 

81  8  128  118  51  129  119  58  8  9 

92  8  129  119  58  130  120  65  9  10 

93  8  130  120  65  131  121  72  10  11 

94  8  131  121  72  132  122  79  11  12 

95  8  132  122  79  133  123  86  12  13 

96  8  133  123  86  134  124  93  13  14 

97  8  134  124  93  135  125  100  14  15 

98  9  136  126  37  137  127  44  6  7 

99  9  137  127  44  138  128  51  7  8 

100  9  138  128  51  139  129  58  8  9 
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4 


101 

9 

139 

129 

58 

140 

130 

65 
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10 

102 

9 

140 

130 

65 

141 

131 

72 

10 

11 

103 

9 

141 
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72 

142 

132 

79 

11 

12 

104 

9 

142 

132 

79 

143 

133 

86 

12 

13 

105 

9 

143 

133 

86 

144 

134 

93 

13 

14 

106 

9 

144 

134 

93 

145 

135 

100 

14 

15 

107 

10 

151 
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37 

152 

137 

44 

6 

7 

108 

10 

152 

137 

44 

153 

138 

51 

7 

8 

109 

10 

153 

138 

51 

154 

139 

58 

8 

9 

110 

10 

154 

139 

58 

155 

140 

65 

9 

10 

111 

10 

155 

140 

65 

156 

141 

72 

10 

11 

112 

10 

156 

141 

72 

157 

142 

79 

11 

12 

113 

10 

157 

142 

79 

158 
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86 

12 

13 

114 

10 

158 

143 

86 

159 

144 

93 

13 

14 

115 

10 

159 

144 

93 

160 

145 

100 

14 

15 

116 

11 

165 

146 

164 

171 

147 

170 

1 

2 

117 

11 

171 

147 

170 

177 

148 

176 

2 

3 

118 

11 

177 

148 

176 

183 

149 

182 

3 

4 

1:9 

11 

183 

149 

182 

189 

150 

188 

4 

5 

1.3 

11 

189 

150 

188 

194 

151 

37 

5 

6 

121 

11 

194 

151 

37 

199 

152 

44 

6 

7 

122 

11 

199 

152 

44 

204 

153 

51 

7 

8 

123 

11 

204 

153 

51 

209 

154 

58 

8 

9 

124 

11 

209 

154 

58 

214 

155 

65 

9 

10 

125 

11 

214 

155 

65 

219 

156 

72 

10 

11 

126 

11 

219 

156 

72 

224 

157 

79 

11 

12 

127 

11 

224 

157 

79 

229 

158 

86 

12 

13 

128 

11 

229 

158 

86 

234 

159 

93 

13 

14 

129 

11 

234 

159 

93 

239 

160 

100 

14 

15 

130 

11 

239 

160 

100 

245 

161 

244 

15 

16 

131 

11 

245 

161 

244 

251 

162 

250 

16 

17 

132 

11 

251 

162 

250 

257 

163 

256 

17 

18 

133 

12 

166 

165 

164 

172 

171 

170 

1 

2 

134 

12 

172 

171 

170 

178 

177 

176 

2 

3 

135 

12 

178 

177 

176 

184 

183 

182 

3 

4 

136 

12 

184 

183 

182 

190 

189 

188 

4 

5 

137 

12 

190 

189 

188 

195 

194 

37 

5 

6 

138 

12 

195 

194 

37 

200 

199 

44 

6 

7 

139 

12 

200 

199 

44 

205 

204 

51 

7 

8 

140 

12 

205 

204 

51 

210 

209 

58 

8 

9 

141 

12 

210 

209 

58 

215 

214 

65 

9 

10 

142 

12 

215 

214 

65 

220 

219 

72 

10 

11 

143 

12 

220 

219 

72 

225 

224 

79 

11 

12 

144 

12 

225 

224 

79 

230 

229 

86 

12 

13 

145 

12 

230 

229 

86 

235 

234 

93 

13 

14 

146 

12 

235 

234 

93 

240 

239 

100 

14 

15 

147 

12 

240 

239 

100 

246 

245 

244 

15 

16 

148 

12 

246 

245 

244 

252 

251 

250 

16 

17 

149 

12 

252 

251 

250 

258 

257 

256 

17 

18 

150 

13 

167 

166 

164 

173 

172 

170 

1 

2 

151 

13 

173 

172 

170 

179 

178 

176 

2 

3 

152 

13 

179 

178 

176 

185 

184 

182 

3 

4 

153 

13 

185 

184 

182 

191 

190 

188 

4 

5 

154 

13 

191 

190 

188 

196 

195 

37 

5 

6 

155 

13 

196 

195 

37 

201 

200 

44 

6 

7 

156 

13 

201 

200 

44 

206 

205 

51 

7 

8 

157 

13 

206 

205 

51 

211 

210 

58 

8 

9 

158 

13 

211 

210 

58 

216 

215 

65 

9 

10 

159 

13 

216 

215 

65 

221 

220 

72 

10 

11 

160 

13 

221 

220 

72 

226 

225 

79 

11 

12 

79 


161 

13 

226 

225 

79 

231 

230 

86 

12 

13 

162  * 

13 

231 

230 

86 

236 

235 

93 

13 

14 

163 

13 

236 

235 

93 

241 

240 

100 

14 

IS 

164 

13 

241 

240 

100 

247 

246 

244 

IS 

16 

165 

13 

247 

246 

244 

253 

252 

250 

16 

17 

166 

13 

253 

252 

250 

259 

258 

256 

17 

18 

167 

14 

168 

167 

164 

174 

173 

170 

1 

2 

168 

14 

174 

173 

170 

180 

179 

176 

2 

3 

169 

14 

180 

179 

176 

186 

185 

182 

3 

4 

170 

14 

186 

185 

182 

192 

191 

188 

4 

5 

171 

14 

192 

191 

188 

197 

196 

37 

5 

6 

172 

14 

197 

196 

37 

202 

201 

44 

6 

7 

173 

14 

202 

201 

44 

207 

206 

51 

7 

8 

174 

14 

207 

206 

51 

212 

211 

58 

8 

9 

175 

14 

212 

211 

58 

217 

216 

65 

9 

10 

176 

14 

217 

216 

65 

222 

221 

72 

10 

11 

177 

14 

222 

221 

72 

227 

226 

79 

11 

12 

178 

14 

227 

226 

79 

232 

231 

86 

12 

13 

179 

14 

232 

231 

86 

237 

236 

93 

13 

14 

180 

14 

237 

236 

93 

242 

241 

100 

14 

15 

181 

14 

242 

241 

100 

248 

247 

244 

15 

16 

182 

14 

248 

247 

244 

254 

253 

250 

16 

17 

183 

14 

254 

253 

250 

260 

259 

256 

17 

18 

184 

15 

169 

168 

164 

175 

174 

170 

1 

2 

185 

15 

175 

174 

170 

181 

180 

176 

2 

3 

186 

15 

181 

180 

176 

187 

186 

182 

3 

4 

187 

15 

187 

186 

182 

193 

192 

188 

4 

5 

188 

15 

193 

192 

188 

198 

197 

37 

5 

6 

189 

15 

198 

197 

37 

203 

202 

44 

6 

7 

190 

15 

203 

202 

44 

208 

207 

51 

7 

8 

191 

15 

208 

207 

51 

213 

212 

58 

8 

9 

192 

15 

213 

212 

58 

218 

217 

65 

9 

10 

193 

15 

218 

217 

65 

223 

222 

72 

10 

11 

194 

15 

223 

222 

72 

228 

227 

79 

11 

12 

195 

15 

228 

227 

79 

233 

232 

86 

12 

13 

196 

15 

233 

232 

86 

238 

237 

93 

13 

14 

197 

15 

238 

237 

93 

243 

242 

100 

14 

15 

198 

15 

243 

242 

100 

249 

248 

244 

15 

16 

199 

15 

249 

248 

244 

255 

254 

250 

16 

17 

200 

15 

255 

254 

250 

261 

260 

256 

17 

18 

200. 

100.  100.  100.  100.  100.  100.  100.  100.  100.  100.  100.  100.  100.  ioo\ft 
100.  100.  f*0 
5  9  5  5 
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APPENDIX  D 
OUTPUT  FILE 

DELTMP  -  200.000000000000000 


I MAT#,  E,  C  MODULI!,  ALFA,  POISSON  AND  HEIGHT  FOLLOW 


1 

0. 10000E+12 

0 . 32051E+11 

0.10000E-03 

0. 30000E+00 

0 . 50000E+00 

2 

0. 10000E+12 

0. 32051E+11 

0. 10000E-03 

0. 30000E+00 

0. 25000E+00 

3 

0.10000E+12 

0.32051E+11 

0. 10000E-03 

0 . 30000E+00 

0. 12500E+00 

4 

0. 10000E+12 

0. 32051E+11 

0. 10000E-03 

0. 30000E+00 

0. 62500E-01 

5 

0.10000E+12 

0. 32051E+11 

0 . 10000E-03 

0 . 30000E+00 

0 . 62500E-01 

6 

0. 10000E+12 

0. 32051E+11 

0.30000E-03 

0. 30000E+00 

0 . 10000E-01 

7 

0. 10000E+12 

0. 32051E+11 

0. 30000E-03 

0. 30000E+00 

0. 20000E-01 

8 

0. 10000E+12 

0. 32051E+11 

0. 30000E-03 

0. 30000E+00 

0 . 40000E-01 

9 

0. 10000E+12 

0. 32051E+11 

0. 30000E-03 

0 . 30000E+00 

0. 20000E-01 

10 

0 . 10000E+12 

0.32051E+11 

0.30000E-03 

0.30000E+00 

0. 10000E-01 

11 

0. 10000E+12 

0. 32051E+11 

0. 10000E-04 

0.30000E+00 

0.62500E-01 

12 

0. 1000QE+12 

0. 32051E+11 

0.10000E-04 

0. 30000E+00 

0.62500E-01 

13 

0.10000E+12 

0. 32051E+11 

0. 10000E-04 

0. 30000E+00 

0.12500E+00 

14 

0 . 10000E+12 

0 . 32051E+11 

0 . 10000E-04 

0.30000E+00 

0. 25000E+00 

15 

0. 10000E+12 

0. 32051E+11 

0 . 10000E-04 

0 . 30000E+00 

0 . 50000E+00 

FORCE 

AND  DISPLACEMENT  VECTORS 

FOLLOW 

1 

0. OOOOOE+OO 

0 . OOOOOE+OO 

2 

0. OOOOOE+OO 

0. OOOOOE+OO 

3 

0.  OOOOOE+OO 

0. OOOOOE+OO 

4 

0. OOOOOE+OO 

0. OOOOOE+OO 

5 

0 . OOOOOE+OO 

0. OOOOOE+OO 
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6 

0 . OOOOOE+OO 

0. OOOOOE+OO 

7 

0.0000024-00 

0. OOOOOE+OO 

8 

0. 00000E+00 

0. 16022E-01 

9 

0. OOOOOE+OO 

-0.51672E-02 

10 

0. OOOOOE+OO 

0. 12042E-01 

11 

0. OOOOOE+OO 

0. 10033E-01 

12 

0. OOOOOE+OO 

0 . 90521E-02 

13 

0. OOOOOE+OO 

0 . 85658E-02 

14 

0. OOOOOE+OO 

0.80331E-02 

15 

0. OOOOOE+OO 

0. 32451E-01 

16 

0. OOOOOE+OO 

-0. 21458E-01 

17 

0. OOOOOE+OO 

0 . 24079E-01 

18 

0. OOOOOE+OO 

0. 19888E-01 

19 

0. OOOOOE+OO 

0. 177712-01 

20 

0. OOOOOE+OO 

0 . 167 16E-01 

21 

0. OOOOOE+OO 

0. 15765E-01 

22 

0. OOOOOE+OO 

0.49296E-01 

23  • 

0. OOOOOE+OO 

-0.50058E-01 

24 

0. OOOOOE+OO 

0 . 36110E-01 

25 

0. OOOOOE+OO 

0.29507E-01 

26 

0. OOOOOE+OO 

0.26255E-01 

27 

0. OOOOOE+OO 

0.24626E-01 

28 

0. OOOOOE+OO 

0.22845E-01 

29 

0. OOOOOE+OO 

0.57798E-01 

30 

0. OOOOOE+OO 

-0.69213E-01 

31 

0. OOOOOE+OO 

0.42198E-01 

32 

0. OOOOOE+OO 

0.34279E-01 

33 

0. OOOOOE+OO 

0. 30238E-01 

34 

0. OOOOOE+OO 

0. 28255E-01 

35 

0. OOOOOE+OO 

0.26546E-01 
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36 

0. OOOOOE+OO 

0.61934E-01 

37 

0. 00000E+00 

-0. 80181E-01 

38 

0. 00000E+00 

0 . 45342E-01 

39 

0. OOOOOE+OO 

0.368UE-01 

40 

0. OOOOOE+OO 

0. 32342E-01 

41 

0. OOOOOE+OO 

0 . 29943E-01 

42 

-0 . 15000E+08 

0 . 27144E-01 

43 

0. OOOOOE+OO 

0.65005E-01 

44 

0. OOOOOE+OO 

-0.88399E-01 

45 

0. OOOOOE+OO 

0. 47763E-01 

46 

0. OOOOOE+OO 

0.38867E-01 

47 

0. OOOOOE+OO 

0.34269E-01 

48 

0. OOOOOE+OO 

0.31904E-01 

49 

0. OOOOOE+OO 

0.29455E-01 

50 

0.  OOOOOE+OO 

0. 68161E-01 

51 

0.  OOOOOE+OO 

-0.97019E-01 

52 

0.  OOOOOE+OO 

0. 50215E-01 

53 

0. OOOOOE+OO 

0.41041E-01 

54 

0. OOOOOE+OO 

0. 36432E-01 

55 

0.  OOOOOE+OO 

0.34170E-01 

56 

0. OOOOOE+OO 

0.32225E-01 

57 

0. OOOOOE+OO 

0.74620E-01 

58 

0. OOOOOE+OO 

-0. 11546E+00 

59 

0.  OOOOOE+OO 

0. 55219E-01 

60 

0.  OOOOOE+OO 

0.45459E-01 

61 

0.  OOOOOE+OO 

0 . 40623E-01 

62 

0. OOOOOE+OO 

0. 38212E-01 

63 

0. OOOOOE+OO 

0.35508E-01 

64 

0. OOOOOE+OO 

0.87846E-01 

65 

0. OOOOOE+OO 

-0.15715E+00 
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66 

0.00000E+00 

0.65382E>01 

67 

0 . 00000E+00 

0.54158E-01 

68 

0.00000E+00 

0.48522E-01 

69 

0 . OOOOOE+OO 

0.45704E-01 

70 

0 . 00000E+00 

0.43179E-01 

71 

0 . 00000E+00 

0. 10123E+00 

72 

0 . 00000E+00 

-0.20527E+00 

73 

0 . 00000E+00 

0.75574E-01 

74 

0, 00000E+00 

0.62744E-01 

75 

0 . 000Q0E+00 

0. 56355E-01 

76 

0 . 00000E+00 

0. 53160E-01 

77 

0 . 00000E+00 

0 . 49681E-01 

78 

O.OOOOOE+OO 

0 . 11433E+00 

79 

0 . OOOOOE+OO 

-0.25976E+00 

80 

O.OOOOOE+OO 

0 . 85778E-01 

81 

O.OOOOOE+OO 

0.71538E-01 

82 

O.OOOOOE+OO 

0. 64345E-01 

83 

O.OOOOOE+OO 

0.60737E-01 

84 

O.OOOOOE+OO 

0.57389E-01 

85 

O.OOOOOE+OO 

0. 12051E+00 

86 

O.OOOOOE+OO 

-0.28934E+00 

87 

O.OOOOOE+OO 

0 . 90866E-01 

88 

O.OOOOOE+OO 

0.76222E-01 

89 

O.OOOOOE+OO 

0.68790E-01 

90 

O.OOOOOE+OO 

0.64988E-01 

91 

O.OOOOOE+OO 

0.60910E-01 

92 

O.OOOOOE+OO 

0. 12339E+00 

93 

O.OOOOOE+OO 

-0 . 30467E+00 

94 

O.OOOOOE+OO 

0.93387E-01 

95 

O.OOOOOE+OO 

0.78682E-01 
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96 

0. 00000E+00 

0.71308E-01 

97 

0 . 00000E+00 

0. 67581E-01 

98 

O.OOOOOE+OO 

0 . 63741E-01 

99 

0 . 25000E+09 

0. 12604E+00 

100 

0.00000E+00 

-0.32032E+00 

101 

0 . 37500E+09 

0.95892E-01 

102 

0.18750E+09 

0 . 81178E-01 

103 

0.93750E+08 

0.73861E-01 

104 

0. 62500E+08 

0.70313E-01 

105 

0.46250E+08 

0.67120E-01 

106 

-0.45000E+08 

0.25831E-01 

107 

O.OOOOOE+OO 

0 . 29799E-01 

108 

O.OOOOOE+OO 

0. 31288E-01 

109 

O.OOOOOE+OO 

0.35675E-01 

110 

O.OOOOOE+OO 

0.42187E-01 

111 

O.OOOOOE+OO 

0.49705E-01 

112 

O.OOOOOE+OO 

0. 56277E-01 

113 

O.OOOOOE+OO 

0. 60826E-01 

114 

O.OOOOOE+OO 

0. 62514E-01 

115 

0 . 45000E+08 

0. 67226E-01 

116 

-0.90000E+08 

0 . 24860E-01 

117 

O.OOOOOE+OO 

0 . 29050E-01 

118 

O.OOOOOE+OO 

0 . 30760E-01 

119 

O.OOOOOE+OO 

0 . 34696E-01 

120 

O.OOOOOE+OO 

0 . 41504E-01 

121 

O.OOOOOE+OO 

0.48460E-01 

122 

O.OOOOOE+OO 

0. 55339E-01 

123 

O.OOOOOE+OO 

0. 59371E-01 

124 

O.OOOOOE+OO 

0.61365E-01 

125 

0.90000E+08 

0 . 65955E-01 
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126 

-0.90000E+08 

0.23563E-01 

127 

0. 00000E+00 

0.27710E-01 

128 

0. 00000E+00 

0.29341E-01 

129 

0 . 00000E+00 

0. 33167E-01 

130 

0 . 00000E+00 

0.39696E-01 

131 

0. 00000E+00 

0. 46422E-01 

132 

0 . 00000E+00 

0.53003E-01 

133 

0 . 00000E+00 

0. 56897E-01 

134 

0 . OOOOOE+OO 

0.58707E-01 

135 

0 . 90000E+08 

0. 63070E-01 

136 

-0.45000E+08 

0. 23197E-01 

137 

0. 00000E+00 

0.27159E-01 

138 

0. OOOOOE+OO 

0. 28411E-01 

139 

0.  OOOOOE+OO 

0.32652E-01 

140 

0.  OOOOOE+OO 

0. 38538E-01 

141 

0.  OOOOOE+OO 

0 . 45660E-01 

142 

0. OOOOOE+OO 

0. 51577E-01 

143 

0.  OOOOOE+OO 

0.55900E-01 

144 

0. OOOOOE+OO 

0.57197E-01 

145 

0 . 45000E+08 

0.61436E-01 

146 

0. OOOOOE+OO 

0. OOOOOE+OO 

147 

0. OOOOOE+OO 

0.71778E-02 

148 

0. OOOOOE+OO 

0.14053E-01 

149 

0. OOOOOE+OO 

0. 20279E-01 

150 

0.  OOOOOE+OO 

0. 23548E-01 

151 

-0. 15000E+08 

0. 23941E-01 

152 

0. OOOOOE+OO 

0.26069E-01 

153 

0. OOOOOE+OO 

0. 28720E-01 

154 

0. OOOOOE+OO 

0.31633E-01 

155 

0. OOOOOE+OO 

0.38716E-01 
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156 

0. 00000E+00 

0.44524E-01 

157 

0 . OOOOUE+OO 

0. 51618E-01 

158 

0 . OOOOOE+OO 

0 . 54641E-01 

159 

0. OOOOOE+OO 

0. 57231E-01 

160 

0. 15000E+08 

0. 59728E-01 

161 

0. OOOOOE+OO 

0 . 59728E-01 

162 

0. OOOOOE+OO 

0.59728E-01 

163 

0.31250E+07 

0 . 59728E-01 

164 

0. OOOOOE+OO 

0. OOOOOE+OO 

165 

0. OOOOOE+OO 

0. OOOOOE+OO 

166 

0. OOOOOE+OO 

0. OOOOOE+OO 

167 

0. OOOOOE+OO 

0. OOOOOE+OO 

168 

0. OOOOOE+OO 

0. OOOOOE+OO 

169 

0. OOOOOE+OO 

0. OOOOOE+OO 

170 

0. OOOOOE+OO 

-0.62365E-02 

171 

0. OOOOOE+OO 

0.66407E-02 

172 

0. OOOOOE+OO 

0.60583E-02 

173 

0. OOOOOE+OO 

0 . 49010E-02 

174 

0. OOOOOE+OO 

0. 26339E-02 

175 

0. OOOOOE+OO 

-0.19392E-02 

176 

0. OOOOOE+OO 

-0.24155E-01 

177 

0. OOOOOE+OO 

0.12867E-01 

178 

0. OOOOOE+OO 

0 . 11783E-01 

179 

0. OOOOOE+OO 

0.96238E-02 

180 

0. OOOOOE+OO 

0. 52627E-02 

181 

0. OOOOOE+OO 

-0.34704E-02 

182 

0. OOOOOE+OO 

-0.52573E-01 

183 

0. OOOOOE+OO 

0. 18851E-01 

184 

0. OOOOOE+OO 

0.17274E-01 

185 

0. OOOOOE+OO 

0 . 14111E-01 
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156 

0 . 00000E+00 

0 . 44524E-01 

157 

0. 00000E+00 

0 . 51618E-01 

158 

0 . 00000E+00 

0.54641E-01 

159 

0 . 00000E+00 

0. 57231E-01 

160 

0. 15000E+08 

0.59728E-01 

161 

0 . 00000E+00 

0 . 59728E-01 

162 

0 . 00000E+00 

0. 59728E-01 

163 

0. 31250E+07 

0. 59728E-01 

164 

0 . 00000E+00 

O.OOOOOE+OO 

165 

0 . 00000E+00 

O.OOOOOE+OO 

166 

0. 00000E+00 

O.OOOOOE+OO 

167 

0. 00000E+00 

O.OOOOOE+OO 

168 

0 . 00000E+00 

O.OOOOOE+OO 

169 

0 . 00000E+00 

O.OOOOOE+OO 

170 

0.00000E+00 

-0.62365E-02 

171 

O.OOOOOE+OO 

0.66407E-02 

172 

0 . 00000E+00 

0.60583E-02 

173 

0.00000E+00 

0 . 49010E-02 

174 

O.OOOOOE+OO 

0 . 26339E-02 

175 

O.OOOOOE+OO 

-0.19392E-02 

176 

O.OOOOOE+OO 

-0.24155E-01 

177 

O.OOOOOE+OO 

0. 12867E-01 

178 

O.OOOOOE+OO 

0. 11783E-01 

179 

O.OOOOOE+OO 

0. 96238E-02 

180 

O.OOOOOE+OO 

0. 52627E-02 

181 

O.OOOOOE+OO 

-0. 34704E-02 

182 

O.OOOOOE+OO 

-0. 52573E-01 

183 

O.OOOOOE+OO 

0. 18851E-01 

184 

O.OOOOOE+OO 

0.17274E-01 

185 

O.OOOOOE+OO 

0. 14111E-01 

L. 
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186 


0 . OOOOOE+OO 


0 . 78855E-02 


187 

0. OOOOOE+OO 

188 

0. OOOOOE+OO 

189 

0. OOOOOE+OO 

190 

0. OOOOOE+OO 

191 

0. OOOOOE+OO 

192 

0. OOOOOE+OO 

193 

0. OOOOOE+OO 

194 

0. OOOOOE+OO 

195 

0.  OOOOOE+OO 

196 

0.  OOOOOE+OO 

197 

0.  OOOOOE+OO 

198 

0.  OOOOOE+OO 

199 

0.  OOOOOE+OO 

200 

0. OOOOOE+OO 

201 

0. OOOOOE+OO 

202 

0. OOOOOE+OO 

203 

0. OOOOOE+OO 

204 

0. OOOOOE+OO 

205 

0. OOOOOE+OO 

206 

0. OOOOOE+OO 

207 

0. OOOOOE+OO 

208 

0. OOOOOE+OO 

209 

0. OOOOOE+OO 

210 

0. OOOOOE+OO 

211 

0. OOOOOE+OO 

212 

0. OOOOOE+OO 

213 

0. OOOOOE+OO 

214 

0. OOOOOE+OO 

215 

0.  OOOOOE+OO 

-0.45865E-02 
-0.70479E-01 
0. 21520E-01 
0 . 19760E-01 
0. 16317E-01 
0 . 92685E-02 
-0.50653E-02 
0.22721E-01 
0.21115E-01 
0. 17566E-01 
0. 10061E-01 
-0.54197E-02 
0.24327E-01 
0.22481E-01 
0 . 18660E-01 
0.10718E-01 
-0.57161E-02 
0. 26229E-01 
0. 24082E-01 
0. 19872E-01 
0.11406E-01 
-0.59275E-02 
0 .  29553E-01 
0.27155E-01 
0.22365E-01 
0.12883E-01 
-0.62042E-02 
0. 35599E-01 
0. 32803E-01 
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216 

O.OOOOOE-fOO 

0.27217E-01 

217 

0. OOOOOE+OO 

0. 15989E-01 

218 

0 . OOOOOE+OO 

-0.64446E-02 

219 

0. OOOOOE+OO 

0 . 4162  IE-01 

220 

0. OOOOOE+OO 

0.38408E-01 

221 

0.  OOOOOE+OO 

0 . 3 1968E-01 

222 

0.  OOOOOE+OO 

0.19138E-01 

223 

0. OOOOOE+OO 

-0.65533E-02 

224 

0. OOOOOE+OO 

0 . 47689E-01 

225 

0. OOOOOE+OO 

0. 44047E-01 

226 

0. OOOOOE+OO 

0.36779E-01 

227 

0. OOOOOE+OO 

0. 22184E-01 

228 

0. OOOOOE+OO 

-0.67637E-02 

229 

0. OOOOOE+OO 

0 . 50965E-01 

230 

0. OOOOOE+OO 

0.47029E-01 

231 

0. OOOOOE+OO 

0. 39072E-01 

232 

0. OOOOOE+OO 

0 . 23484E-01 

233 

0. OOOOOE+OO 

-0. 68603E-02 

234 

0. OOOOOE+OO 

0.52884E-01 

235 

0. OOOOOE+OO 

0 . 48468E-01 

236 

0. OOOOOE+OO 

0 . 39986E-01 

237 

0. OOOOOE+OO 

0. 23982E-01 

238 

0. OOOOOE+OO 

-0. 68389E-02 

239 

0. OOOOOE+OO 

0. 54036E-01 

240 

0. OOOOOE+OO 

0 . 49242E-01 

241 

0. OOOOOE+OO 

0.40488E-01 

242 

0. OOOOOE+OO 

0 . 24321E-01 

243 

0. OOOOOE+OO 

-0.66809E-02 

244 

0. OOOOOE+OO 

-0. 32032E+00 

245 

0. OOOOOE+OO 

0. 54036E-01 
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246 

0. 00000E+00 

0.49242E-01 

247 

0 . 00000E+00 

0.40488E-01 

248 

0 . 00000E+00 

0 . 2432  IE-01 

249 

0. 00000E+00 

-0 . 66809E-02 

250 

0 . 00000E+00 

-0.32032E+00 

251 

0 . 00000E+00 

0 . 54036E-01 

252 

0. 00000E+00 

0. 49242E-01 

253 

0 . 00000E+00 

0 . 40488E-01 

254 

0. 00000E+00 

0 . 24321E-01 

255 

0 . 00000E+00 

-0.66809E-02 

256 

0 . OOOOOE+OO 

-0. 32032E+00 

257 

0. 62500E+07 

0.54036E-01 

258 

0 . 93750E+07 

0.49242E-01 

259 

0. 18750E+08 

0. 40488E-01 

260 

0. 37500E+08 

0. 24321E-01 

261 

0.25000E+08 

-0.66809E-02 

STRESSES  FOLLOW 

SGMABOT 

SGMATOP  TAUELEM 

1 

-0.96832E+08 

0. 20163E+09  0 . 33346E+07 

2 

-0.97221E+08 

0 . 23223E+09  0.42668E+07 

3 

-0. 97662E+08 

0.26336E+09  0.34617E+07 

4 

-0.86943E+08 

0.27526E+09  0.17059E+07 

5 

-0. 56664 E+08 

0. 24077E+09  -0.22757E+08 

6 

-0.31570E+08 

0.22852E+09  0.30736E+08 

7 

~0. 19144E+08 

0 . 26242E+09  0.22740E+08 

8 

0.80498E+06 

0 . 29174E+09  0.14828E+08 

9 

0. 16305E+08 

0.32262E+09  0.54632E+07 

10 

0.19130E+08 

0. 33807E+09  0.15863E+06 
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11 

0 . 204 16E+08 

0 . 30984E+09 

-0.93561E+07 

12 

0. 17685E+08 

0. 23665E+09 

-0.30987E+08 

13 

0.85421E+07 

0. 15424E+09 

-0 . 54002E+08 

14 

0. 19906E+07 

0 . 58011E+08 

-0.78675E+08 

15 

“0 . 24752E+09 

-0 . 96832E+08 

0 . 45844E+07 

16 

-0 . 26088E+09 

*0. 97221E+08 

0. 58611E+07 

17 

-0.27857E+09 

-0. 97662E+08 

0. 44885E+07 

18 

-0.28418E+09 

-0.86943E+08 

0. 99650E+07 

19 

-0. 24033E+09 

-0.56664E+08 

-0.13647E+06 

20 

-0.17785E+09 

-0 . 31570E+08 

0 . 63424E+08 

21 

-0 . 13044E+09 

-0 . 19144E+08 

0. 53303E+08 

22 

“0. 11640E+09 

0.80498E+06 

0.31608E+08 

23 

-0. 13006E+09 

0. 16305E+08 

0 . 883 19E+07 

24 

-0. 14138£«-09 

0. 19130E+08 

-0 . 16488E+06 

25 

“0. 12063E+09 

0.20416E+08 

“0. 11367E+08 

26 

-0.63188E+08 

0.17685E+08 

-0.44518E+08 

27 

-0.16015E+08 

0 . 85421E+07 

-0 . 84443E+08 

28 

-0.17644E+07 

0. 19906E+07 

-0.12074E+09 

29 

-0.32109E+09 

-0.247E2E+09 

0. 15406E+07 

30 

-0.34612E+09 

-0. 26088E+09 

0 . 55932E+07 

31 

-0 . 36365E+09 

-0. 27857E+09 

0. 86522E+06 

32 

-0 . 40251E+09 

-0 . 28418E+09 

0. 14019E+08 

33 

-0.36879E+09 

-0 . 24033E+09 

0.36426E+08 

34 

-0.22946E+09 

-0.17785E+09 

0. 10874E+09 

35 

-0 . 13485E+09 

“0. 13044E+09 

0 . 75305E+08 

36 

-0.16177E+09 

-0 . 11640E+09 

0 . 28705E+08 

37 

-0. 21003E+09 

-0 . 13006E+09 

0.61505E+07 

38 

-0. 21678E+09 

-0.14138E+09 

-0.35792E+06 

39 

-0. 20095E+09 

-0.12063E+09 

-0. 52206E+07 

40 

-0. 11105E+09 

-0.63188E+08 

-0. 21053E+08 
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41 

0 .  70878E+07 

-0.16015E+08 

-0.67557E+08 

42 

0 . 2 13  04E+08 

-0.17644E+07 

-0. 12319E+09 

43 

-0.35756E+09 

-0 . 32109E+09 

0. 48163E+06 

44 

-0.38872E+09 

-0.34612E+09 

0. 34010E+07 

45 

-0.40678E+09 

-0.36365E+09 

0 . 57895E+06 

46 

-0.45566E+09 

-0 . 40251E+09 

0.54280E+07 

47 

-0.49365E+09 

-0 . 36879E+09 

0 . 69360E+08 

48 

-0.21563E+09 

-0. 22946E+09 

0 . 16808E+09 

49 

-0.93571E+08 

-0.13485E+09 

0 . 81349E+08 

50 

-0. 19161E+09 

-0. 16177E+09 

0 . 16037E+08 

51 

-0.25080E+09 

-0.21003E+09 

0 . 46495E+07 

52 

-0.2 54 3 3 E+09 

-0.21678E+09 

-0 . 24469E+06 

53 

-0.24238E+09 

-0. 20095E+09 

-0.22418E+07 

54 

-0 . 14975E+09 

-0. 11105E+09 

0.41097E+07 

55 

0. 37217E+08 

0.70878E+07 

-0. 35371E+08 

56 

0.92679E+08 

0.21304E+08 

-0. 14127E+09 

57 

-0.39752E+09 

-0. 35756E+09 

0. 12391E+08 

58 

-0 . 42008E+09 

-0.38872E+09 

-0 . 11199E+08 

59 

-0. 46904E+09 

-0.40678E+09 

0. 13011E+08 

60 

-0 . 44486E+09 

-0. 45566E+09 

-0 . 26006E+08 

61 

-0 . 82061E+09 

-0 . 49365E+09 

0 . 10136E+09 

62 

-0.75409E+08 

-0 . 21563E+09 

0 . 29180E+09 

63 

0. 10776E+09 

-0.93571E+08 

0. 21460E+08 

64 

-0. 34326E+09 

-0. 19161E+09 

0.97283E+07 

65 

-0. 23293E+09 

-0 . 25080E+09 

0.41638E+07 

66 

-0 . 34978E+09 

-0.25433E+09 

-0.26277E+07 

67 

-0 . 22926E+09 

-0. 24238E+09 

0. 39777E+07 

68 

-0 . 29578E+09 

-0 . 14975E+09 

0. 80823E+07 

69 

0.13236E+09 

0. 37217E+08 

0 . 64530E+08 

70 

0. 35154E+09 

0.92679E+08 

-0.20324E+09 
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71 

“0. 14128E+10 

~0 . 20754E+10 

0. 50013E+09 

72 

-0 . 24042E+10 

“0. 18922E+10 

-0.15450E+09 

73 

~0 . 21226E+10 

-0.23433E+10 

0. 51903E+08 

74 

-0 . 23488E+10 

-0.22329E+10 

-0.13578E+08 

75 

-0 . 22482E+10 

-0.23498E+10 

0.97283E+07 

76 

-0 . 23428E+10 

-0. 22293E+10 

-0 . 26196E+07 

77 

-0 . 20902E+10 

-0.22958E+10 

0. 20390E+08 

78 

-0. 23248E+10 

-0.18676E+10 

0. 13S14E+09 

79 

-0. 11151E+10 

-0 . 16485E+10 

-0.21072E+09 

80 

-0. 13241E+10 

-0. 14128E+10 

0.32412E+09 

81 

-0 . 23160E+10 

-0. 24042E+10 

-0 . 81814E+08 

82 

-0.22127E+10 

-0. 21226E+10 

0 . 25195E+08 

83 

-0 . 23192E+10 

-0.23488E+10 

-0.50890E+07 

84 

-0.23045E+10 

-0. 22482E+10 

0 . 29591E+07 

85 

~0 . 23 120E+10 

-0 . 23428E+10 

0 . 26092E+07 

86 

-0.21937E+10 

-0 . 20902E+10 

0 . 21170E+08 

87 

-0 . 22022E+10 

-0. 23248E+10 

0. 12030E+09 

88 

-0 . 11639E+10 

-0.U151E+10 

-0. 68495E+08 

89 

-0.13412E+10 

-0. 13241E+10 

0 . 30606E+07 

90 

-0.23475E+10 

-0. 23160E+10 

0. 42052E+06 

91 

-0.22348E+10 

-0.22127E+10 

-0 . 89273E+06 

92 

“0. 23470E+10 

-0. 23192E+10 

0. 68946E+06 

93 

-0.23275E+10 

-0 . 23045E+10 

-0.12984E+07 

94 

-0. 23418E+10 

-0.23120E+10 

0.57419E+07 

95 

-0.22213E+10 

-0.21937E+10 

0 . 30953E+08 

96 

-0 . 22758E+10 

-0.22022E+10 

0.90213E+08 

97 

-0. 12551E+10 

-0. 11639E+10 

0 . 21441E+09 

98 

-0.14152E+10 

-0.13412E+10 

-0. 31989E+09 

99 

-0. 24992E+10 

-0.23475E+10 

0.80553E+08 

100 

-0.21520E+10 

-0 . 22348E+10 

-0.24553E+08 
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101 

-0. 24114E+10 

-0.23470E+10 

102 

-0.22878E+10 

-0.23275E+10 

103 

-0. 24083E+10 

-0.23418E+10 

104 

-0.21353E+10 

-0.22213E+10 

105 

“0. 24811E+10 

-0.22758E+10 

106 

-0.13045E+10 

-0 . 12551E+10 

107 

-0.21488E+10 

-0. 14152E+10 

108 

-0.19394E+10 

-0.24992E+10 

109 

-0 . 24 174E+10 

-0. 21520E+10 

110 

-0 . 22918E+10 

-0.24114E+10 

111 

-0. 24192E+10 

-0. 22878E+10 

112 

-0 . 22906E+10 

-0.24083E+10 

113 

-0. 23955E+10 

-0 . 21353E+10 

114 

-0 . 19640E+10 

-0 . 24811E+10 

115 

-0.20010E+10 

-0. 13045E+10 

116 

0.39805E+09 

0.43834E+09 

117 

0. 36699E+09 

0. 41563E+09 

118 

0. 34874E+09 

0. 36697E+09 

119 

0 . 30036E+09 

0.39032E+09 

120 

0 . 26053E+09 

0.17822E+08 

121 

0. 54227E+09 

0.75123E+09 

122 

0.66086E+09 

0.96063E+09 

123 

0. 56485E+09 

0.48262E+09 

124 

0 . 50460E+09 

0. 60825E+09 

125 

0.50214E+09 

0. 48083E+09 

126 

0 . 50685E+09 

0. 60941E+09 

127 

0 . 55520E+09 

0. 50449E+09 

128 

0 . 66764E+09 

0.93602E+09 

129 

0.36045E+09 

0 . 89902E+09 

130 

59468. 

-0. 12464E+06 

0 . 48791E+07 
“0. 36850E+07 
0 . 66497E+07 
0.45573E+08 
0 . 42590E+08 
0. 51209E+09 
-0.49752E+09 
0. 14742E+09 
-0. 45783E+08 
0 . 99768E+07 
-0.66053E+07 
0 . 74915E+07 
0 . 55718E+08 
-0.10628E+07 
0. 67659E+09 
-0. 12190E+08 
0. 10995E+08 
-0.12801E+08 
0 . 25549E+08 
-0. 10015E+09 
-0.29445E+09 
-0. 19738E+08 
-0.10095E+08 
-0.39153E+07 
0. 14889E+07 
0. 53087E+07 
0.53776E+08 
0 . 91013E+08 
0.56745E+09 
0 . 79085E+09 
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131 

35680. 

-74781. 

132 

11893. 

-24927. 

133 

0.35437E+09 

0 . 39B05E+09 

134 

0. 32933E+09 

0 . 36699E+09 

135 

0. 31182E+09 

0 . 34874E+09 

136 

0 . 27291E+09 

0 . 30036E+09 

137 

0 . 30660E+09 

0 . 26053E+09 

138 

0.44624E+09 

0 . 54227E+09 

139 

0. 54043E+09 

0 . 66086E+09 

140 

0. 51459E+09 

0. 56485E+09 

141 

0 . 46481E+09 

0. 50460E+09 

142 

0. 46054E+09 

0. 50214E+09 

143 

0 .  46389E+09 

0 . 50685E+09 

144 

0 . 49645E+09 

0 . 55520E+09 

145 

0.47539E+09 

0 . 66764E+09 

146 

0 . 20980E+09 

0. 36045E+09 

147 

-2635.9 

59468. 

148 

-1521.7 

35680. 

149 

-507.45 

11893. 

150 

0.26757E+09 

0.35437E+09 

151 

0. 25423E+09 

0 . 32933E+09 

152 

0 . 23653E+09 

0.31182E+09 

153 

0 . 23091E+09 

0 . 27291E+09 

154 

0 . 27481E+09 

0 . 30660E+09 

155 

0 . 33732E+09 

0.44624E+09 

156 

0 . 38S09E+09 

0 . 54043E+09 

157 

O.39860E+09 

0.51459E+09 

158 

0. 38514E+09 

0.46481E+09 

159 

0.37509E+09 

0.46054E+09 

0 . 79085E+09 
0.79085E+09 
-0. 57935E+06 
-0.33031E+07 
-0.67715E+06 
-0.52327E+07 
-0 . 69720E+08 
-0.16835E+09 
-0 . 81001E+08 
-0. 16459E+06 
-0.41742E+07 
-0.12006E+07 
0 .  11140E+08 
0. 47156E+08 
0. 17594E+09 
0.35491E+09 
0 . 32965E+09 
0 . 32965E+09 
0. 32965E+09 
-0.15447E+07 
-0. 55885E+07 
-0.87117E+06 
-0.14003E+08 
-0.36549E+08 
-0.10886E+09 
-0.75559E+08 
-0. 28515E+08 
-0.61843E+07 
-0 . X9486E+06 


160 


0. 38117E+09 


0.46389E+09 


0. 10904E+08 


161 

0. 35850E+09 

0. 49645E+09 

162 

0 . 26551E+09 

0.47539E+09 

163 

0.10077E+09 

0. 20980E+09 

164 

3473.6 

-2535.9 

165 

2084.3 

-1521.7 

166 

694.89 

-507.45 

167 

0.97544E+08 

0 . 26757E+09 

168 

0 . 97157E+08 

0 . 25421E+09 

169 

0.96713E+08 

0 . 236S3E+09 

170 

0.1074SE+09 

0.23091E+09 

171 

0. 13776E+09 

0. 27481E+09 

172 

0.16289E+09 

0 . 33732E+09 

173 

0. 17526E+09 

0. 38509E+09 

174 

0. 19540E+09 

0 . 39860E+09 

175 

0.21053E+09 

0. 38514E+09 

176 

0.21496E+09 

0.37509E+09 

177 

0 . 20458E+09 

0. 38117E+09 

178 

0 . 15989E+09 

0. 35850E+09 

179 

0. 99276E+08 

0.26551E+09 

180 

0.35805E+08 

0. 10077E+09 

181 

1175.8 

3473.6 

182 

705.45 

2084.3 

183 

235.10 

694.89 

184 

-0.24544E+09 

0.97544E+08 

185 

-0. 21484E+09 

0.97157E+08 

186 

“0. 18370E+09 

0. 96713E+08 

187 

-0. 17182E+09 

0. 10745E+09 

188 

-0.20633E+09 

0.13776E+09 

189 

-0. 21854E+09 

0 . 16289E+09 

190 

-0. 18458E+09 

0. 17526E+09 

0. 55894E+08 
0 . 14 181E+09 
0 . 203 19E+09 
0.  11629E+09 
0. 11629E+09 
0.  11629E+09 
-0.45917E+07 
-0. 58539E+07 
-0.44958E+07 
-0.99519E+07 
61981. 

-0.63509E+08 
-0. 53328E+08 
-0. 31672E+08 
-0.87828E+07 
81373. 

0. 11422E+08 
0.38791E+08 
0. 59307E+08 
0. 55522E+08 
-0 . 55906E+08 
-0. 5S906E+08 
-0.55906E+08 
-0.33428E+07 
-0 . 42584E+07 
-0.34705E+07 
-0.16839E+07 
0. 22719E+08 
-0. 30764E+08 
-0. 22765E+08 
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191 

-0.15534E+09 

0.19540E+09 

192 

-0.12404E+09 

0 . 21053E+09 

193 

-0.11087E+09 

0. 21496E+09 

194 

~0 . 12 104E+09 

0 . 20458E+09 

195 

-0.11932E+09 

0 . 15989E+09 

196 

-0.91445E+08 

0. 99276E+08 

197 

-0.3677SE+08 

0.35805E+08 

198 

-2705.9 

1175.8 

199 

-1623.6 

705.45 

200 

-541.10 

235.10 

UPPER 

PEEL 

STRESS 

m 

-0. 13210E+09 

UPPER 

PEEL 

STRESS 

m 

-0. 53388E+08 

UPPER 

PEEL 

STRESS 

m 

-0 . 19007E+08 

UPPER 

PEEL 

STRESS 

m 

-0 . 64454E+07 

UPPER 

PEEL 

STRESS 

m 

-0.78777E+07 

UPPER 

PEEL 

STRESS 

m 

-0. 27239E+08 

UPPER 

PEEL 

STRESS 

m 

-0.69977E+08 

UPPER 

PEEL 

STRESS 

m 

-0 . 20688E+09 

LOWER 

PEEL 

STRESS 

— 

-0 . 13335E+09 

LOWER 

PEEL 

STRESS 

m 

-0.53074E+08 

LOWER 

PEEL 

STRESS 

m 

-0.19075E+08 

LOWER 

PEEL 

STRESS 

m 

-0 . 64680E+07 

LOWER 

PEEL 

STRESS 

m 

-0.73789E+07 

LOWER 

PEEL 

STRESS 

m 

-0 . 23411E+08 

LOWER 

PEEL 

STRESS 

m 

-0 . 56982E+08 

LOWER 

PEEL 

STRESS 

m 

-0 . 15128E+09 

-0 . 1482 1E+08 
-0. 55446E+0? 
0. 41476E+06 
0 . 47076E+07 
0.43322E+07 
-0.54208E+07 
-0.25159E+08 
-0. 14 118E+09 
-0 . 14 116E+09 
-0. 14118E+09 
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